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ABSTRACT 


A  Laser  Doppler  Velocimeter  (LDV)  was  used  to  perform  a 
propeller  wake  survey  for  DTNSRDC  propellers  4381  and  4383. 

Three  components  of  the  velocity  were  measured  and  from  this, 
the  velocity  field  and  the  vorticity  field  were  constructed. 
Through  a  coordinate  transformation,  the  vorticity  field  was 
separated  into  two  parts  i.e.,  trailing  vorticity  and  boundary 
layer  vorticity  in  the  viscous  wake. 

The  data  was  used  for  the  prediction  of  propeller  blade 
profile  drag  (viscous  sectional  drag) .  The  idea  of  hypothetical 
flow  by  Betz  was  extended  to  propeller  flow  and  a  formula  for 
the  profile  drag  was  derived  in  terms  of  velocities,  which  is 
suitable  for  LDV  application.  The  radial  distribution  of 
profile  drag  was  computed  using  the  velocity  data.  It  was 
shown  that  the  value  of  Cd  is  close  to  that  of  2-dimensional 
section  at  mid-radius  range,  but  becomes  substantially  lower 
at  inner  and  outer  range,  but  becomes  substantially  lower 
at  inner  and  outer  radii  and  that  it  goes  up  again  near  the 
tip  and  the  hub.  This  trend,  is  consistent  with  the  behavior 
of  profile  drag  of  a  finite  span  wing. 
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1.  Introduction 

Marine  screw  propeller  is  relatively  new  in  the  long 
history  of  ships  which  is  comparable  to  the  history  of 
mankind  itself.  It  was  only  a  hundred  years  ago  when  the 
propeller  took  the  present  style[l].  Since  then,  through 
various  evolutionary  stages,  it  has  been  used  almost 
exclusively  as  a  mean  of  ship  propulsion. 

During  last  two  decades,  rapid  growth  of  computer 
technology  has  made  it  possible  for  theoretical  and 
numerical  methods  to  analyze  various  problems  in  science  and 
engineering.  The  field  of  marine  propeller  hydrodynamics  is 
not  an  exception.  More  than  a  few  computer  programs  have 
been  or  are  being  developed  using  numerical  lifting-line  or 
lifting-surface  theory. 

At  MIT,  there  has  been  a  continuous  effort  for  over  20 
years  by  Kerwin  and  his  colleagues  to  develop  con^puter 
programs  for  the  design  of  propeller  blades  and  for  the 
prediction  of  propeller  performance,  based  on  numerical 
lifting  surface  theory[2].  The  effort  has  produced  some  of 
the  best  propeller  programs  around,  such  as  PUF-2[3], 
PUF-3C4],  PBD-10C5]  and  FPV-10[6]. 

At  almost  the  same  period  of  time.  Laser  Doppler 
Anemometry  was  introduced  as  a  powerful  tool  in  experimental 
fluid  dynamics.  In  certain  areas  of  experiments,  it  enabled 
us  to  read  flow  velocity  directly  with  a  very  good  accuracy 
and  spatial  resolution  without  disturbing  the  flow.  At  the 
Department  of  Ocean  Engineerig,  MIT,  Laser  Doppler 
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Velocimeter  (LDV)  has  been  in  operation  since  1977.  Min[7] 
was  the  first  to  set  up  the  Laser  apparatus  and  bring  it 
into  work.  He  showed  various  posibilities  of  the  flow 
measurement  around  a  propeller,  thus  proving  the 
powerfulness  of  LDV.  After  having  successfully  shown  its 
capability,  the  LDV  is  now  being  used  to  back  up  those 
computer  programs  as  well  as  in  various  experiments. 

This  thesis  is  intended  to  extend  the  work  done  by  Min. 
Among  several  things  to  be  done,  the  author  has  been 
particularly  interested  in  the  prediction  of  the  profile 
drag  (viscous  sectional  drag)  of  a  propeller  blade,  using 
the  velocity  data  obtained  by  LDV,  which  has  never  been 
tried  so  far.* 

The  current  versions  of  the  programs  use  a  constant 
value  for  the  coefficient  of  the  profile  drag  (Cd)  over  the 
radius  of  the  propeller  blade.  By  feeding  the  actually 
measured  values  of  the  drag  distribution,  it  is  expected 
that  the  accuracy  of  the  programs  will  improve,  especially 
of  the  design  program  (PBD-10).  One  of  the  objectives  of 
this  paper,  therefore,  is  to  supply  such  information  using 
the  experimental  data  obtained  by  the  LDV. 

In  this  paper,  the  method  of  Betz[8]  for  the  prediction 
of  profile  drag  which  uses  the  concept  of  hypothetical  flow 


*) There  have  been  some  examples  of  LDV  application  to  lift 
measurement  such  as  Sayre[9]  or  OrlofftlO].  But  as  far  as 
the  author  knows,  there  is  no  example  of  profile  drag 
measurement  by  LDV. 
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will  be  generalized  first  and  a  formula  will  be  given  in 
terms  of  velocities  as  opposed  to  in  terms  of  pressures 
which  was  given  by  Betz.  A  method  to  obtain  the 
hypothetical  flow  from  measured  velocity  data  will  also  be 
discussed. 

The  measured  velocity  field,  and  vorticity  distribution 
which  is  derived  from  the  velocity  data  will  be  presented 
also,  and  the  flow  field  in  the  propeller  wake  will  be 
discussed.  The  propellers  tested  are  DTNSRDC  research 
propeller  s  4381  (no  blade  skew)  and  4383  (with  72  deg  skew) . 
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2.  Experimental  Facility 

MIT  Water  Tunnel  The  experiments  were  done  in  the  Variable 
Pressure  Water  Tunnel  of  the  Marine  Hydrodynamics 
Laboratory,  Department  of  Ocean  Engineering,  MIT. 

As  Min[7]  mentions,  this  water  tunnel  is  ideal  for 
Laser  measurements.  The  tunnel  is  a  closed  circulating 
channel,  having  a  square  test  section  with  transparent 
plastic (plexiglass)  viewing  windows  on  four  sides.  The  size 
of  the  test  section  is  20"  x  20"  and  the  viewing  window  is 
44"  x  16"  x  2"  thick.  The  tunnel  operates  in  the  speed 
range  of  0  to  30  ft/sec. 

LDV  system  The  LDV  system  is  manufactured  by  Thermo  Systems 
Inc.  with  Spectra  Physics  15mW  He-Ne  Laser  and  has  Dual 
Beam  Forward  Scatter  Mode.  The  LDV  mount  has  a  traverse 
system  which  has  3-degrees  of  freedom.  Traversing  is  done 
by  turning  the  handles  by  hand,  and  the  position  of  the 
Laser  beam  crossing  point  is  obtained  by  reading  the  scales 
attached  to  the  mount.  This  traverse  gear  is  rather  simple. 
However,  careful  reading  of  the  scales  can  give  as  much  as 
0.1mm  resolution. 

Data  Acquisition  Data  acquisition  and  reduction  is  done  by 
a  mini-conputer  (Digital  Equipment  MINC-11)  through  an  A/D 
converter.  Data  sampling  is  possible  at  a  rate  as  high  as 
20kHz.  If  the  propeller  rotates  at  900rpm,  for  instance, 
this  gives  a  resolution  of  0.27  degree  in  angle.  The 
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characteristics  of  the  MIT  Water  Tunnel  LDV  and  Data 
acquisition  systems  are  summarized  in  Table  2-1.  The 
experimental  set-up  is  illustrated  in  Figure  2-1  and  the 
diagram  of  data  acquisition  system  is  illustrated  in  Figure 
2-2. 


Propellers  The  propellers  tested  are  DTNSRDC  4381  and  4383 
research  propellers  of  12  inch  diameter,  out  of  a  series  of 
5  propellers  with  systematically  changing  skew  and 
rake[ll,12].  Propeller  4381  is  the  parent  propeller  of  the 
series  and  has  no  blade  skew.  Propeller  4383  has  72  degree 
blade  skew  and  also  skew  induced  rake.  The  principal 
dimensions  are  given  in  Table  2-2.  They  were  designed  to 
give  the  same  performance  characteristics  (KT,Kg,n)  at  the 
design  condition  of  J=.889.  Figures  2-3  through  2-7  show 
the  result  of  open  water  tests  done  in  the  MIT  Water  Tunnel 
recently.  The  propeller  efficiencies  at  design  J  are 
similar  to  each  other  except  for  Prop. 4381  which  has  lower 
efficiency  than  the  others.  This  is  due  to  the  fact  that, 
as  the  parent  propeller,  it  has  been  used  more  than  any 
others  in  the  series  and  has  been  wearing  and  getting 
roughened  leading  edges. 

Wall  Ef feet  For  experiments  in  the  water  tunnel,  there  is 
always  a  propblem  to  be  kept  in  mind.  That  is  the  tunnel 
wall  effect  or  so-called  blockage  effect.  The  presence  of 
the  walls  give  a  restriction  to  the  flow.  The  streamlines 
on  the  wall  have  to  be  parallel  to  it,  since  no  flow  into 


the  wall  is  possible.  This  distorts  the  flow  from  what  it 
would  be  in  an  unbounded  fluid  flow.  The  significance  of 
this  effect  varies  depending  on  the  kind  and  size  of  the 
test  model  one  uses,  as  well  as  the  flow  conditions.  For 
the  propeller  4381  operating  at  the  design  condition,  the 
magnitude  of  the  wall  effect  was  checked  using  FPV-10.  By 
replacing  the  walls  with  the  images  of  the  propeller,  it  is 
possible  roughly  to  know  the  wall  effect.  At  J=.889  and  at 
0.7R,  for  example,  this  effect  is  about  a  few  per  cent  at 
most.  Hence  it  can  be  said  that  the  wall  effect  is 
negligible  in  this  particular  case.  In  the  computation  of 
profile  drag,  the  flow  is  assumed  to  be  unbounded. 
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3. Experiments  and  Results 


3 . 1  Field  Point  Velocity  and  Vorticity 

Coordinate  Systems  Four  coordinate  systems  are  used  : 

i)  Cartesian  coordinates  o-xyz 

where  the  x-axis  is  along  the  propeller  shaft  center  line, 
pointing  upstream.  The  y-axis  is  set  horizontally  toward 
starboard  and  the  z-axis  is  vertically  downward.  The  y-z 
plane  corresponds  to  the  center  plane  of  propeller  rotation. 

ii)  Cylindrical  coordinate  system  o-xre 
where  r  is  radially  outward,  0  goes  clockwise  when  looked 
from  downstream,  and  the  r-e  plane  is  identical  to  y-z 
plane.  0=0  corresponds  to  the  y-axis. 

There  is  a  relation  between  the  two  systems,  which  is 


z~  rsine  J 

The  above  two  are  space-fixed  (or  tunnel-fixed)  frames  of 
reference.  The  following  are  blade  attached  coordinate 
systems  : 

iii)  Cartesian  coordinates  o-xy*z* 

iv)  Cylindrical  coordinates  o-xra ' 

Each  of  the  two  moving  systems  is  defined  to  be  in  the  same 
way  as  i)  or  ii)  at  time  0. 

The  relations  between  the  coordinate  systems  are  : 


0' 

=  0  -  at 

] 

■ 

y' 
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where  n  is  angular  velocity  of  the  propeller  rotation, 
which  is  constant  in  this  paper. 

The  unit  vectors  are  also  defined  as  shown  in  the 
figure  below. 


; 

Figure  3-1  Coordinate  systems 


Construction  of  Velocity  Field  The  Laser  system  in  the  MIT 
Water  Tunnel  is  a  one  component  Laser  which  reads  one 
component  of  the  velocity  at  a  time.  Therefore,  to  get 
three  components,  three  identical  running  conditions  have  to 
be  repeated.  The  axial  component  is  taken  essentially  at 
any  location  as  long  as  the  beam  crossing  point  is  visible 
from  the  photo- multiplier.  But  it  is  convenient  if  the 
measurement  point  is  on  the  z-axis  where  the  radial 
component  can  be  also  measured.  The  radial  component  is 
taken  as  the  vertical  component  on  z-axis,  and  the 
tangential  conponent  can  be  measured  as  the  vertical 
velocity  component  along  the  y-axis.  The  velocity  field  can 
be  constructed  from  the  3  velocity  components  thus  obtained. 
This  technique  is  only  valid  when  the  flow  around  the 
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propeller  is  steady,  in  other  words,  constant  propeller 
rotation  in  a  uniform  flow  or  in  an  axisymmetrical  flow. 
The  phase  lag  between  the  velocity  components  measured  at 
different  places  will  be  taken  care  of  by  shifting  the  phase 
of  the  data  in  the  computer. 

The  propellers  were  run  at  900rpm  most  of  the  time, 
which  gives  blade  Reynolds  No.  of  0.5  xl0‘  to  1.5  xl0‘.  A 
few  data  were  taken  at  600  and  1200rpm  also.  One  limitation 
is  that  we  can  not  get  broad  range  of  Reynolds  No.  in  MIT 
Water  Tunnel  if  we  are  to  run  the  propeller  at  one  fixed  J. 

The  Laser  signal  usually  contains  noise  from  various 
sources  as  mentioned  by  Min[7].  Removal  of  noise  from  the 
real  signal  is  necessary.  One  possibility  is  averaging 
(ensemble  average) .  If  the  flow  is  laminar,  this  technique 
gives  good  results.  If  the  flow  contains  turbulence, 
however,  this  will  eliminate  the  tu-fulencr  component 
together  with  the  noise.  Taking  rms  curbulence  may  be 
satisfactory,  but  discrimination  of  noise  from  turbulence 
remains  a  problem.  Although  it  is  said  that  turbulence 
measurement  is  possible  by  LDV[13],  the  author  is  not  quite 
convinced  of  that,  at  least  with  the  propeller  flow.  When 
calculating  drag  from  wake  survey  data,  velocity  comes  in 
quadratic  form  in  the  momentum  flux.  If  there  is  10% 
turbulence  in  the  flow,  for  example,  this  becomes  a  1% 
fluctuation  in  momentum.  Fortunately,  in  most  of  the 
conditions  discussed  in  this  paper,  the  turbulence  level  in 
the  propeller  wake  is  considered  to  be  relatively  small. 
Hence,  turbulence  is  neglected  and  only  mean  velocities  are 
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measured. 

Figure  3-5  shows  an  example  of  ensemble  average  of  the 
velocity  for  360  degrees. 

Figures  3-6  a,b,c  to  3-24  a,b,c  show  the  field  point 
velocities  taken  at  various  radii  on  the  x=-0.333R 
(downstream)  plane  for  Prop. 4381  and  on  the  x=-0.  488R  plane 
for  Prop. 4383.  Each  of  these  positions  (in  x)  is  almost  the 
closest  possible  to  each  propeller  in  the  downstream 
direction. 

3  components  'taken 


I 


Figure  3-2  Restriction  to  the  LDV  measurement 

In  table  3-1  is  shown  the  distance  of  each  measurement 
point  from  the  trailing  edge,  along  the  helical  line 
constructed  with  uniform  inflow  and  propeller  rotation. 

Data  sampling  frequency  was  either  10kHz  or  20kHz  and 
ensemble  average  of  1000  or  2000  was  taken  for  each 
component.  Only  one  blade  passage  is  given  here. 

In  these  figures,  sub  number  a  is  for  x-component 
(axial),  b  is  for  r-component  (radial)  and  c  is  for 
e-component  (tangential).  Note  that  the  frame  of  reference 
is  fixed  in  space  so  that  e -component  does  not  include 
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virtual  velocity  due  to  propeller  rotation  which  will  be 
observed  from  propeller-fixed  frame  of  reference.  The  angle 
in  the  abscissa  corresponds  to  -9  ' ,  since  the  measurement 
point  moves  in  the  negative  e' -direction  against  the 
propeller.  Smaller  angles  correspond  to  the  pressure  side 
of  the  blade. 

Note  :  The  propeller  rpm  for  these  figures  is  900rpm  except 
for  Figures  3-12,  3-22  (600rpm)  and  Figure  3-13  (1200rpm). 

Vorticity  Field  Vorticity  of  the  flow  is  given  by 

HlP  =  y  x  D  (3-3) 


Since  the  velocities  are  easily  obtained,  the  computation  of 
vorticity  can  be  done  through  the  formula  (3-3)  numerically. 
If  we  discretize  Eq.(3-3),  we  have 


(3-4) 


Where  i  is  the  subscript  for  x,  j  is  for  r  and  k  is  for  e . 

To  do  this  confutation,  four  more  measurement  stations 
were  set  around  the  field  point  of  interest  and  velocities 


were  measured. 
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Figure  3-3  Measurement  stations  for  Vorticity 


Figure  3-25  shows  an  example  of  thus  computed  vorticities  at 

the  measurement  point  x*=-0.333R,  r=0.  7R  for  Propeller  4381. 

For  the  convenience  of  coordinate  transformation  which  will 

be  discussed  later  in  drag  calculation,  the  frame  of 

reference  is  set  to  the  propeller  blade  and  the  computation 

of  vorticity  was  done  in  xre*  coordinates.  Therefore,  twice 

the  propeller  rotation,  2a,  is  introduced  into  x-component. 

To  distinguish  the  vorticity  in  blade  attached  coordinates 

from  that  in  fixed  frame  of  reference,  the  notations  co^,  ujy 
* 

and  ojg-  are  used.  If  0J«  ,  60*.  and  coe  represent  the 
vorticity  in  space  fixed  coordinates,  the  relation  between 
the  two  vorticities  is  s 

2^  j  <-or  ~  cOr^  a '  — 


3 . 2  Discussion  on  Velocity  and  Vorticity  Field 
Velocity  Field  The  three  components  of  the  velocity 
and  vorticity  tell  us  various  things.  First  of  all,  if  we 
look  at  the  radial  component  of  velocity,  we  find  that  the 
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velocity  jump  across  the  trailing  vortex  sheet  is  there. 
(In  the  real  flow,  of  course,  the  sheet  is  not  really  a 
sheet  of  zero  thickness,  tout  has  a  finite  thickness.  The 
velocity  jump  has  a  width  for  transition  as  a  consequence. 
However,  the  word  'sheet'  will  be  used  in  this  paper  for 
convenience. )  Although  this  component  is  not  all  for  the 
flow  due  to  trailing  vortex,  it  is  considered  to  be  very 
close  to  it,  especially  for  the  no  skewed  propeller 
(Prop. 4381).  Vortex  core  type  flow  is  well  represented 
there  as  we  expect  (Figure3-llb  for  example) . 

This  velocity  jump  switches  its  sign  (the  shape  becomes 
the  other  way)  as  we  go  along  r-direction.  For  Prop. 4381, 
2Jr  at  r/R=0.6  and  above  are  different  from  that  at  r/R=0.5 
and  below.  The  implication  is  that  the  slope  of  the  bound 
circulation  changes  its  sign  and  consequently  the  sign  of 
trailing  vortex  sheet  changes  at  a  point  between  the  two 
corresponding  radii.  This  occurs  on  Prop. 4383  between 
r/R=0.6  and  0.7  (Figures  3-20b  and  3-21b) . 

As  the  measurement  point  moves  radially  outward  toward 
the  tip,  the  influence  of  tip  vortex  becomes  pronounced. 
Since  vorticity  in  the  tip  vortex  is  much  stronger  than  the 
trailing  vortex  sheet,  the  velocity  jump  becomes  much 
larger( Fig.3-15b  and  3-24b)  .  zr%  and  also  pick  up  the 
influence  of  tip  vortex.  The  large  hump  in  Vx  (Fig.3-15a 
and  3-24a)  is  the  flow  due  to  the  tip  vortex.  Since  this 
measurement  point  is  still  inside  the  slip  stream  (0.91R  for 
Prop. 4381,  0.95R  for  Prop. 4383  at  these  measurement 
locations),  it  appears  as  strong  negative  x-flow  (in  other 
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words,  the  same  direction  as  freestream) .  At  a  point 
outside  the  propeller  wake,  the  tip  vortex  gives  positive 
x-flow  as  shown  in  Figure  3-16a.  When  this  comes  to  2J» ,  it 
is  a  little  confusing,  because  it  appears  as  a  dent  as  shown 
in  Figs.  3-15c  and  24c.  The  first  (i.e.  at  smaller  angle) 
of  the  two  dents  in  either  figure  is  considered  to  be  the 
wake  which  is  due  to  the  boundary  layer  and  the  second  is 
the  flow  due  to  tip  vortex.  It  is  interesting  to  see  that 
this  conqponent  (e-conponent)  hardly  picks  up  the  influence 
of  tip  vortex  once  outside  the  wake  (Fig.3-16c). 

The  wake  of  the  blade  (which  comes  off  the  boundary 
layer  on  the  blade)  appears  as  velocity  defect  (or  dent) 
typically  in  x-  or  8 -component  at  most  of  the  radii  of 
measurement.  (Let  us  call  this  'viscous  wake'.)  However,  by 
carefully  examining  the  radial  velocity  component,  we  find 
that  it  is  also  there.  The  little  hump  on  top  of  the 
velocity  peak  on  pressure  side  of  Ur  in  Figure  3-llb  for 
example,  is  considered  to  be  a  projection  of  viscous  wake. 
Across  the  sign  change  of  trailing  vortex,  it  then  appears 
on  the  suction  side  as  in  Figure  3-9b.  Since  the  flow  field 
is  under  the  strong  influence  of  slipstream  contraction,  we 
can  not  say  that  the  hump  always  has  outward  velocity. 
However,  it  consistently  shows  a  tendency  to  go  outward 
compared  with  the  other  part.  This  can  be  explained  by 
centrifugal  force.  Due  to  the  no-slip  condition,  the  fluid 
in  the  boundary  layer  on  the  blade  tends  to  move  together 
with  the  blade.  When  it  leaves  the  blade  to  form  the  wake, 
it  tries  to  flow  along  the  tangent  line  to  that  radius. 
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giving  radial  component  to  the  velocity.  It  is  also 
consistent  with  the  idea  of  centrifugal  force  that  the 
radial  speed  of  the  little  hump  increases  toward  outer 
radii.  The  magnitude  of  outward  velocity  may  not  change  in 
proportion  to  r.  It  is  expected  that  at  some  radii  on  the 
blade,  there  is  a  transition  from  laminar  to  turbulent  flow. 
Once  it  is  turbulent,  the  flow  should  have  less  influence 
from  the  blade[14],  changing  the  magnitude  of  outward  flow 
in  the  wake. 

In  Figures  3-26  and  3-27  are  shown  the  projection  of 
the  velocity  defect,  at  each  radius  of  measurement,  to  the 
direction  of  the  flow  in  the  viscous  wake  U -direction, 
say).  The  coordinate  system  here  is  blade  attached.  Hence, 
the  velocity  is  larger  at  outer  radii.  Note  that  the 
abscissa  is  blade  passage  (in  non-dimensional 
distance,  r/Rde 1 )  and  not  in  angle,  consequently  giving 
different  lengths  to  one  blade  passage  between  different 
radii.  The  velocity  defect  is  deeper  at  mid-radius  range 
(in  the  vicinity  of  r/R=0.6),  where  the  flow  is  considered 
to  be  reasonably  2-dimensional.  Near  the  blade  tip,  or  near 
the  propeller  hub,  the  velocity  defect  becomes  shallower. 
This  is  particularly  so  for  Prop. 4383.  One  of  the  possible 
reasons  for  this  is  that  there  is  a  nontrivial  3-dimensional 
effect.  At  outer  radii,  relatively  large  centrifugal  force 
and  roll-up  of  the  tip  vortex  may  be  causing  cross  flow. 
Near  the  propeller  hub,  where  the  distance  between  the 
blades  becomes  smaller  and  where  blades  go  into  the  hub, 
strong  interactions  between  them  possibly  exist,  creating 
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another  3-dimensional  effect. 

One  interesting  thing  is  that  the  two  sheets  -  trailing 
vortex  and  viscous  wake  -  are  not  in  the  same  location,  as 
might  be  noticed  when  looking  at  the  little  hump  in  The 
hump  is  clearly  off  the  center  of  velocity  jump.  This  is 
not  surprising  at  all  since  the  origins  of  the  two  wakes  are 
different.  The  trailing  vortex  is  created  by  the  presence 
of  tip  vortex  while  viscous  wake  comes  from  the  boundary 
layers  on  the  blade.  Figures  3-28  and  3-29  illustrate  the 
geometry  of  the  two  sheets  and  also  in-plane  velocity  of 
each  point  on  the  sheets  for  the  two  propellers  observed  on 
the  measurment  planes  (x=-0.333R  for  Prop. 4381  and  x=-0.488R 
for  Prop. 4383) .  These  graphs  were  obtained  by  finding  the 
angular  location  of  the  viscous  wake  and  the  center  of  the 
velocity  jump  using  the  graphs  of  and  .  The  graphs 
clearly  show  the  different  locations  and  motions  of  the  two 
sheets.  The  sheets  for  Prop. 4383  have  larger  curvature  than 
Prop. 4381,  implying  the  effect  of  large  skew  of  the  blade. 
At  outer  radii,  the  trailing  vortex  sheet  is  fairly  far  away 
from  the  viscous  wake,  toward  suction  side.  At  the  inner 
radii,  the  two  sheets  almost  line  up( .4R  and  below) . 
Correspondingly,  the  viscous  wake  does  not  show  itself  in 
Vr  • 

Field  point  velocities  were  also  measured  at  x=-1.0R 
and  -2. OR  at  r=0.7R  for  Prop. 4381.  They  are  shown  in 
Figures  3-30  and  3-31.  Two  blade  passages  are  shown  in 
these  figures.  At  x»-1.0R,  the  velocity  defect  and  velocity 
jump  are  still  there  although  becoming  shallower  and  smaller 


-31- 


because  of  diffusion.  But  at  x»*-2.0R  they  seem  to  have 
disappeared.  It  would  be  interesting  to  examine  what  is 
going  on  between  these  two  points  by  taking  finer  spacing  in 
x  and  getting  more  data. 

Comparison  with  FPV-10  output  As  a  byproduct  of  PUF-2,  a 
field  point  velocity  program  (FPV-10[6])  was  written  and  has 
been  in  use  to  predict  field  point  velocities  around  the 
propeller.  It  is  now  possible  by  using  the  LDV  to  map  the 
velocity  field  and  back  up  the  program. 

At  some  measurement  points,  velocity  components  were 
compared  with  FPV-10  output.  Since  the  program  uses  a 
potential  flow  model,  we  should  not  expect  the  velocity 
defect  in  the  wake.  First,  comparisons  were  made  in  the 
potential  flow  region  i.e.  outside  the  wake.  They  are 
shown  in  Figures  3-32  and  3-33.  There  is  some  D.C.  offset 
between  the  two,  especially  in  the  x-components .  The 
amplitudes  in  the  radial  components  are  a  little  discrepant. 
Tangential  components  are  close  to  each  other.  Phases  in 
each  velocity  component  are  in  good  agreement.  A 
modification  to  this  program  is  being  made  as  PSFFPV  by 
Greeley[15].  An  output  of  this  program  for  the  same  field 
point  is  shown  in  Figure  3-34.  Although  there  is  D.C. 
offset,  the  amplitude  of  each  component  is  closer  to  the 
measured  velocity. 

Once  the  field  point  is  in  the  wake,  there  is  some 
difference  (Figs.  3-35  and  3-36).  While  the  measured 
velocity  shows  a  velocity  defect  in  the  wake  of  the  blade, 
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the  confuted  velocity  does  not.  Since  both  velocities 
(measured  and  computed)  must  satisfy  the  continuity 
equation,  discrepancy  between  the  two  is  larger  to 
compensate  the  velocity  defect  in  the  wake. 

Figure  3-37  for  the  conputed  velocity  (Prop. 4383 
x=-.47R,  r=0.924R)  shows  the  influence  of  tip  vortex  fairly 
well  as  coipared  with  Figure  3-38,  the  measured.  The  peak 
in  Vn.  (influence  by  tip  vortex)  in  the  computed  velocity 
field  is  sharper  than  that  measured.  The  difference  is 
presumably  due  to  the  presence  of  the  boundary  layer  in  the 
real  flow.  It  is  considered  that  the  boundary  layer  smears 
the  tip  vortex  to  some  extent. 

Vorticity  Field  The  author  tried  to  separate  the  vorticity 
in  the  wake  of  the  propeller  into  two,  corresponding  to  the 
two  different  types  of  wake,  namely  trailing  vortex  wake  and 
viscous  wake.  This  idea  is  also  used  by  Tan[16]  in  the 
analysis  of  the  wake  of  gas  turbine  rotor.  The  vorticity 
distributions  obtained  by  Eq.(3-4),  one  example  being  shown 
in  Figure  3-25  were  transformed  to  a  new  coordinate  system 
o'-5nC«  where  c-axis  is  parallel  to  the  flow  at  the  center 
of  trailing  vortex  wake.  The  second  axis  n  was  made  to  lie 
on  the  trailing  vortex  sheet.  Consequently,  the  third  axis 
C  points  normal  to  the  trailing  vortex  sheet.  The 
transformed  vorticities  are  shown  in  Figures  3-39  through 
3-56. 

Since  5-axis  is  set  parallel  to  the  trailing  vortex, 
aJ>£  should  represent  trailing  vorticity.  This  seems  to  be 
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almost  successful  at  most  radii.  However,  as  discussed  in 
the  previous  section,  there  is  a  component  of  viscous  wake 
projected  to  r-direction  and  thus  to  the  direction  of  the 
trailing  vorticity.  Therefore  it  is  not  completely  possible 
to  separate  trailing  vorticity  from  the  others  only  by 
coordinate  transformation.  This  projection  appears  as  a 
small  bump  at  the  foot  of  large  hump  of  trailing  vorticity 
as  shown  in  Figure  3-44  for  instance.  Except  for  that,  cOf 
is  a  fairly  good  even  function  of  6. 

n-coraponent  shows  typical  vorticity  distribution  in  the 
viscous  wake.  Since  the  wake  comes  from  the  boundary  layers 
on  both  sides  of  the  blade,  this  vortex  component  has  two 
vortices  with  opposite  signs.  Hence  the  vorticity  is  an  odd 
function  of  9 .  The  third  component  (  U> ^  )  shows  very  little 
vorticity  as  is  expected  through  the  coordinate 
transformation.  Due  to  the  curvature  of  the  wake  sheets  and 
also  due  to  the  fact  that  the  two  wake  (trailing  vortex  and 
viscous)  sheets  are  a  little  apart,  this  coitponent  does  not 
become  completely  zero,  or  flat  except  at  a  few  radii. 


Figure  3-4  Coordinate  Transformation 


The  vorticity  plot  becomes  more  noisy  at  inner  radii 
toward  the  propeller  hub.  This  is  presumed  ma’nly  due  to 
the  finer  spacing  of  data  sampling.  Since  the  sampling 
frequency  was  kept  constant,  in  other  words  angle  interval 
was  constant,  it  gave  smaller  spacing  between  the  data 
points  at  smaller  radii.  The  smaller  the  interval  between 
the  points  is,  the  larger  the  error  in  discretized  vorticity 
becomes.  But  U)$  comes  out  as  a  nice  even  function  of  e 
there . 

At  r=0.875R  for  both  propellers,  the  vorticities  show 
fairly  complicated  aspect,  especially  in  n  and  ; 
components.  This  is  considered  to  be  due  to  the  strong 
influence  of  the  tip  vortex.  In  certain  region  on  the  blade 
near  the  tip,  a  part  of  the  boundary  layer  is  pulled 
laterally  (radially)  toward  the  tip  vortex  and  enters  the 
swirling  flow  just  outside  the  tip  vortex  core.  And  it 
makes  'roll  cake'  type  structure  together  with  the  vortex 
core.  This  phenomenon  was  observed  more  distinctly  in  the 
wake  of  a  hydrofoil.  (See  Section  4.3) 

Aerofoil  theory  tells  us  that  the  strength  of  the 
trailing  vorticity  is  proportional  to  the  slope  of  bound 
circulation  on  the  blade.  The  area  under  cO$  has  a 
dimension  of  circulation/length  and  corresponds  to  the  slope 
of  circulation.  At  each  radius,  the  integration  was  done 
and  thus  obtained  slopes  are  plotted  in  Figure  3-57  for  the 
both  propellers.  From  the  analogy  to  an  airplane  wing,  one 
might  expect  for  these  curves  to  decrease  monotonically  from 
inner  radii  to  outer  radii.  However,  in  the  case  of  these 
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propellers,  there  is  a  decrease  of  the  value  toward  the  hub. 
This  trend  is  considered  to  be  due  to  the  presence  of  the 
hub. 

In  the  same  figure  are  also  plotted  the  slopes  (dr/dr) 
obtained  by  PUF-2  program.  They  are  shown  by  a  solid  line 


and  a  dashed  line.  Agreement  between  the  two  (measured 

and 

computed) 

is  good  at 

mid-radius 

range.  But 

near 

the 

propeller 

hub,  PUF-2 

overpredicts 

the  slope. 

This 

is 

partially  because  PUF-2  neglects  the  hub.  If  the  hub  gives 
some  2-dimensional  effect  as  does  the  water  tunnel  wall  on  a 
2-D  hydrofoil,  the  slope  should  tend  to  approach  zero.  The 
measured  data  implies  that.  Near  the  tip  on  the  other  hand, 
PUF-2  is  underpredicting  the  slope. 

Error  in  Vorticity  The  vorticity  field  was  obtained  by 
numerically  differentiating  the  velocity.  Differentiation 
usually  enlarges  the  error  (or  noise)  which  is  contained  in 
the  original  signal.  Aside  the  discretization  error,  the 
followings  are  the  possible  sources  of  error  : 

Error  in  positioning 
Error  in  velocity  reading. 

Let  us  discuss  how  these  come  into  the  error  of  the 
vorticity.  Take  the  radial  component  as  an  example. 
Discretized  vorticity  is  written  as 


Now  suppose  that  each  element  contains  an  error,  like 


A  &  -  A©  C  I  ■+  £„e 
A  %  =  Aa°(  |  t  E**) 

where  superscript  0  represents  true  value  and  £  with 
subscript  is  the  error.  Expand  the  velocity  components  in 
terms  of  errors.  We  have 


Mg,.,  =  vlk„  +  e*vlt„  +  “-e 


3  O 


=  ^U-  +  C  ulk.,  +  A© 


where  the  second  term  on  the  right  of  each  equation  is  the 
reading  error  and  the  third  is  the  positioning  error.  Other 
quantities  can  be  expressed  in  the  same  way.  The  vorticity 
then  becomes 


CO.  =.  - - -  |?;0L  -  7r°,  +<Cr  7  r®  _  C"  7  r°  .  r-  *  Afl 
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If  we  make  rms  error,  we  have 


yms  error  = 


|oo^.|  1  N  Y*B 


I  dVxk 
y  a© 


S'X 


What  we  need  at  this  point  is  an  estimate  of  error  in 
each  element.  The  error  is  due  to  the  positioning. 
That  is  about  10%  with  the  current  traverse  system.  £A0  is 
mainly  due  to  the  error  in  the  computer  clock  which  is 
negligibly  small.  The  error  in  velocity  reading  can  be 
small  also. 

A  rough  estimation  showed  an  error  on  the  order  of'  10% 
which  mainly  comes  from  the  positioning.  The  current 
positioning  system  is  rather  crude  as  mentioned  in  Chapter 
2.  A  modification  to  the  traverse  and  positioning  system  is 
being  planned.  It  is  expected  that  the  error  is  reduced  by 
this. 

Fluctuation  of  the  vorticity  itself  (so-called  vortex 
wandering[17 ] )  may  also  cause  an  error.  This  might  be  a 
matter  of  reproducibility  rather  than  error,  however. 
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4. Profile  Drag  of  Propeller  31ade  Section 

The  velocity  and  vorticity  data  obtained  from  the 
experiments  and  shown  in  Chapter  3  will  be  used  in  the 
computation  of  propeller  blade  profile  drag  (viscous 
sectional  drag) .  To  do  sc,  a  formula  for  the  drag 
coefficient  has  to  be  developed. 

A  handy  formula  for  the  profile  drag  of  an  aerofoil  was 
first  derived  by  Betz[8].  To  get  the  force  on  a  body  in  a 
flow,  momentum  theorem  has  to  be  used.  To  do  so,  one  has  to 
measure  velocity  and  pressure  all  over  the  control  surface 
which  surrounds  the  body.  Noticing  the  fact  that  the  flow 
outside  the  aerofoil  wake  is  practically  inviscid 
irrotational,  Betz  introduced  a  hypothetical  flow  which  is 
inviscid  and  irrotational  and  which  has  the  same  pressure  as 
the  real  flow  all  over  the  control  surface  and  the  same 
velocity  on  the  control  surface  except  in  the  wake.  Thus  he 
eliminated  the  need  of  measuring  the  flow  except  in  the 
wake.  His  formula  is  expressed  in  terms  of  pressures.  This 
was  adequate  because  what  one  measured  then  was  pressure 
(total  and  static) . 

With  the  rapid  development  of  Laser  Doppler  Anemometry, 
direct  velocity  measurement  is  now  possible.  Unlike  the 
pressure  probe  such  as  Pitot  tube,  LDV  does  not  disturb  the 
flow.  High  frequency  response  of  LDV  is  also  suitable  for 
measurements  of  rapidly  changing  flow  such  as  propeller 
flow. 
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The  author's  intention  here  is  to  derive  a  formula 
which  gives  profile  drag  of  propeller  blade,  and  which  is 
written  in  terms  of  velocities,  thus  suitable  for  LDV 
application. 

In  this  chapter,  a  formula  to  give  viscous  drag  on  a 
solid  body  in  a  uniform  flow  is  derived  by  extending  and 
generalizing  Betz's  method.  Then  it  is  applied  to 
hydrofoils  and  propellers. 

4.1  Derivation  of  the  Formula 

Moment un  Theorem  The  force  and  moment  which  act  on  a  solid 
body  in  a  flow  are  obtained  by  considering  a  momentum 
balance.  Suppose  there  is  a  solid  body  moving  in  an 
incompressible  fluid  flow.  Place  a  control  volume  V' 
(material  volume)  around  the  body  at  a  particular  time,  t 
say. 


Figure  4-1  Formulation  of  the  problem 
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Newton*  s  second  law  of  motion  gives 

+  ffif  JdTT-  F  (4-1. ) 

T'  4*  T' 

P^dV=-f^AP‘ylld^'l'©rxf3^V'W-B  (4-2) 

■V'  S'  ▼' 


density  of  the  fluid 
velocity  vector 
surface  of  control  volume  V' 
surface  of  the  body 
static  pressure 

normal  unit  vector  on  S*  and  SB 
pointing  outward  the  control  volume 
acceleration  of  gravity 
position  vector 

the  force  which  the  body  receives  from  the  flow 
the  moment  "  " 

The  fluid  receives  a  reaction  from  the  body  as  -  and 
-&v  Here,  shearing  stresses  on  the  control  surface  s'  due 
to  viscosity  of  the  fluid  were  neglected  because  they  are 
considered  to  be  small  even  in  the  wake  of  the  body  at  most 
of  the  conditions  in  which  we  are  interested  in. 

Since  the  flow  measurement  is  done  at  points  fixed  in 
space,  it  is  advantageous  to  use  space  fixed  control  volume. 
If  use  is  made  of  transport  theorem[18]  on  Equations  (4-1) 
and  (4-2),  the  material  volume  is  replaced  by  space  fixed 
volume.  And  if  we  rearrange  them  to  get  expressions  for 


where 

P 

V 

J$' 

P 

to 


r  : 
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force  and  moment,  we  have 


F0=  -%(?VV*+Wd&*%?%Av-\\\^dv-\\<?VVM 

£  ▼  r  St 

and 

-  JJj  ^  P|f^  V  -  ^  W> *  P  VUnd$ 

*1 

where  V  and  S  are  fixed  in  space. 

For  a  solid  body,  surface  integrals  over  SB  in  both 
equations  are  zero  since  there  is  no  net  momentum  transfer 
across  S3.  Thus,  we  obtain 

(4-3) 

J  TT  "V 

{4-4) 

$  v  -r 

These  two  equations  are  the  ones  which  basically  give 
the  force  and  moment  from  measured  velocity  and  pressure. 
However,  these  formulas  are  not  so  practical.  The  reason  is 
obvious.  Velocity  and  pressure  have  to  be  evaluated  all 
over  the  control  surface.  And  if  the  flow  is  unsteady,  the 
acceleration  term  has  to  be  evaluated  everywhere  in  the 
volume.  At  this  point,  it  is  necessary  to  introduce 
hypothetical  flow. 

Introduction  of  Hypoth et ical  flow  The  presence  of  the  solid 
body  introduces  two  different  types  of  flows.  Downstream 
the  body,  there  is  a  strip  of  flow  in  which  the  influence  of 
the  body  and  of  the  viscosity  of  the  fluid  are  not 


Vorticity  which  is 


negligible.  This  is  called  wake, 
generated  in  the  boundary  layer  on  the  body  is  convected 
downstream  and  forms  the  wake.  Therefore,  the  wake  is 
highly  rotational.  If  one  dimension  of  the  body 
perpendicular  to  the  flow  is  small,  like  a  hydrofoil,  the 
wake  is  usually  confined  to  a  narrow  region.  Outside  the 
wake,  the  flow  remains  almost  the  same  as  is  in  the  incoming 
flow.  For  example,  if  the  inflow  is  uniform,  then  the 
outside  flow  is  inviscid  and  irrotational.  There, 
Bernoulli's  theorem  holds. 

Betz  extended  inviscid  irrotational  flow  into  the  wake 
of  the  aerofoil  as  hypothetical  flow  in  order  to  use 
Bernoulli's  theorem.  However,  in  a  3-D  flow  such  as  that 
around  a  loaded  airplane  wing  or  propeller,  Bernoulli's 
theorem  may  not  hold  because  the  flow  in  the  wake  is 
rotational  due  to  trailing  vortex.  Therefore,  the  author 
does  not  claim  that  the  hypothetical  flow  is  irrotational. 
Even  so,  it  is  still  advantageous  to  introduce  inviscid 
hypothetical  flow. 

Let  velocity  Tp  and  static  pressure  be  those  in 

the  hypothetical  flow.  This  hypothetical  flow  is  to  be 
created  by  placing  distribution  of  singularities  such  as 
sources  and  vortices  in  the  incoming  flow.  As  previously 
mentioned,  the  idea  is  to  eliminate  the  need  of  measuring 
velocity  and  pressure  all  over  the  control  surface. 
Therefore,  it  is  necessary  that  the  hypothetical  flow  has 
the  same  velocity  and  pressure  distribution  on  S  as  the  real 
flow  at  least  in  the  inviscid  flow  region. 
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The  requirements  for  the  hypothetical  flow  are: 

The  hypothetical  flow  is  inviscid  flow. 

On  the  control  surface  S,  pressure  p*  is  everywhere 
identical  to  pressure  p  of  the  real  flow  and  velocity 
sf'  is  identical  to  D  everywhere  except  in  the  wake. 


The  two  flows  become  different  in  the  wake  because  one 
is  viscous  flow  and  the  other  is  inviscid.  This  means  that 
there  is  extra  flux  in  the  hypothetical  flow. 


Figure  4-2  Hypothetical  flow 


The  force  and  moment  on  the  singularities  due  to  the 
flow  are  obtained  by  again  applying  momentum  balance 
argument, 


—  -  f’  <«-s> 

V*'  $’  v*' 

+  $  AfxPJchr  -JtlT  (4-6) 

sv*'  S'  v*' 


where  ff*~  and  ifcIL*  are  the  force  and  moment  on  the  source  and 
vortex  distribution.  V’  represents  the  control  volume  for 
the  hypothetical  flow,  which  is  a. little  different  from  V' 

x' 

in  that  V  does  not  have  the  solid  body  in  it  but  the 

*  / 

singularities  instead.  The  outer  surface  of  V  is  taken  to 
be  identical  to  S’  here. 

It  is  necessary  to  bring  Equations  (4-5)  and  (4-6)  into 
the  same  form  as  Equations  (4-3)  and  (4-4).  Since  the 
velocity  has  singularities,  care  has  to  be  taken  in 
manipulating  the  integrals  which  involve  velocity. 

The  hypothetical  flow  consists  of  3  different  types  of 
flows;  free  stream  (LU,  flow  from  the  sources  TU*  and  flow 
due  to  vortices  <■ 

7U*-  (LU  +  V>/  +  V*  (4“7) 

For  the  vortex  part  of  the  flow,  vortex  core  type 
representation  is  desired  rather  than  vortex  line.  As  will 
be  mentioned  later,  the  hypothetical  flow  is  asked  to 
simulate  lifting  problem  of  the  body.  This  means  that  the 
vorticity  pattern  in  the  hypothetical  flow  which  is  related 
to  the  lifting  problem  has  to  be  the  same  as  that  in  the 
real  flow.  One  possible  place  where  this  vorticity  pattern 
is  seen  is  the  wake  on  the  control  surface  S,  where  the 
trailing  vortex  passes  through  toward  downstream.  As  Figure 
3-llb  shows,  it  is  typically  vortex  core  type  flow,  which 
makes  sense  since  there  is  no  vortex  line  (i.e.  zero 
diameter)  in  the  real  flow.  £j^then  becomes  continuous  and 
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is  not  singular.  The  only  singularity  comes  from  2J>S  •  Hence 
special  care  is  needed  when  integrating  ^  .  The  integral 
is  of  Cauchy  type  and  is  determined  by  first  excluding  an 
infinitesimal  volunief usually  sphere)  which  contains  the 
source,  and  taking  the  limit  of  zero  radius,  like 


U \v\fo  =  (4-8) 

V*' 

Here,  V£  represents  infinitesimal  volume  (sphere) 

surrounding  the  source  at  ( **>/&,  Z* )  • 

Using  this,  the  left  side  of  (4-5)  becomes 

?u>*dv  =  +  *£  )dv~ 

V#'  "V* ' 

=  w„+a>i-)<<v  +  ?M\  v?dlr 

V*  V*' 

= f  +  pjj(uviC)iC«is 

y*  >s 

a.*,*  ri’9C^e 

=  p®^.(U.t®0<<T+  pj(cy<.+jj*)u-*d5 
▼  *  £ 

+ P&M  n*$  +P*r®(fe;  §&«#} 

where  (p^  is  the  velocity  potential  for  the  source 
distribution.  The  velocity  potential  does  not  exist  for 
Here,  the  transport  theorem  and  Gauss1  theorem  were 
used  in  the  manipulation. 

The  third  term  on  the  right  reduces  to 

&  &  s 


I 
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This  relation  is  derived  if  we  use  space-fixed  surface  S  in 
(4-5)  and  conpare  the  result  with  that  obtained  by  using 
material  surface. 

Then, 

T*'  -r*  $ 


Finally,  Equation  (4-5)  reduces  to 


d'Xs^s^s 


(4-9  ) 


The  moment  (Equation  (4-6))  can  be  handled  in  a  similar 

way . 

)dr 

v*'  W 

=  fAt  $  ir*  (U«>^  )<*v  +$  rwtfdY] 

T*'  -V*' 


But 


r>  v<?5/  =  - 


can  be  shown  by  vector  calculus.  And  if  we  use  the  relation 


?&§  $(rm)d$  -  ?\\  -f^Kor*  **)<*$  +  pjj  ( f*  &>/)  v*  a$ 

S'  s  s 

which  is  obtained  again  by  comparing  the  formulas  for 
angular  moment un  with  2  different  types  of  control  surfaces. 
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one  moving  with  the  fluid,  the  other  fixed  in  space.  Then, 


we  get 


a*  %  ( or*  r*f)dv  =  ?%  nr*^((u^iV)4Tr  *p$$ir*aviV)i#*£ 

tt*'  t*  5 

+  P5\^^ITx«HS  +  IT*  2D*  V>*d$ 

&  5 

*t&$)\cMgV»n)<iS )  dtsd%dis 

w.  «•"’& 


We  end  up  with 


tt*  =  -ftr*  ?«&£■*•?«  MS  *  P  85  r*gdv-  pjk&Cu+i&rfir 

/S  V+  v* 

-eS#W)4$  -  (4-10) 


w* 


Decomposit ion  of  Force  and  Moment  The  forces  and  moments 
are  decomposed  as  follows. 

For  the  real  flow,  we  have 


fb  - 


(4-11) 


tv+  Iti.o+  Mtwor 


(4-12) 


where  is  the  lifting  force  produced  by  inflow  +  induced 

velocity,  is  the  drag  force  and  is  the  buoyancy 

force.  These  notions  also  apply  to  moment.  Note  that 
'lift*  here  means  the  force  which  circulation  distribution 
around  the  body  produces  against  inflow  and  induced 
velocity.  Consequently,  this  includes  so-called  induced 
drag.  This  convention  is  the  same  in  the  hypothetical  flow. 
The  force  w  Lch  the  source  distribution  receives  is  not 


called  ’lift*  in  this  paper  even  though  the  force  may  happen 
to  be  parallel  to  the  lift. 


It  is  clear  that 

Uw  =  -^r»Pfdv 

For  the  hypothetical  flow  we  can  write 

IF*  -  fff  +  Fo 

M*  -  mu  +  M.0 


(4-13) 

(4-14) 

(4-15) 

(4-16) 


*  * 

where  F.  is  the  lifting  force  and  TFo  is  the  drag.  There  is 
no  buoyancy  force  in  the  hypothetical  flow. 

Since  we  art  interested  in  the  drag,  it  would  be  most 
convenient  if  if:  -  if*.  and  (MU  =  ItHgv.  were  true.  One 
necessary  condition  for  this  is  that  the  vorticity 
distribution  in  the  hypothetical  flow  is  equal  to  that  which 
contributes  to  lift  in  the  real  flow.  This  will  be  discused 
later  with  each  specific  case. 

Subtracting  Equations  (4-9)  and  (4-10)  from  (4-3)  and 
(4-4)  respectively  and  taking  into  account  the  relations 
(4-11)  -  (4-16),  we  have 

**♦«*-£- IF? 

Sw  v  v* 

+  (4-n) 

S  w. 


% 
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“V +  -  n:  -  n' 

=  -?^Df*(TL)L ?\\W*£:(V*>+V*r)ciV- 

+  f))| £(r<«14S  4>*ar*tn)d$)ch,cl$sclzs  (4-18) 

where  Sw  represents  the  part  of  the  control  surface  S  which 
is  in  the  wake  of  the  body. 


The  Force  and  Moment 


in  Hypothetical  Flow  The  net 


flux  from  the  sources  is  not  zero  because  velocity  TO  is 
different  from  TO  in  the  wake.  The  real  flow  satisfies  the 
continuity  equation.  The  net  flux  produces  a  non-zero  force 
on  the  sources  (Lagally  force) .  If  we  write  the  flux  Q  in 
terms  of  source  strength,  q(  ),  we  have 


£==$$$  &('Xs,3s,Zs)<=<'X»^2, 
(or,  =  SMCV.^Vi  ) 


(4-19) 


If  the  sources  are  moving  and  also  changing  their 
intensity  in  time,  q  becomes  a  function  of  time.  But  in 
this  paper,  let  us  just  assume  no  time  variation  of 
intensity  and  no  relative  motion  between  the  sources  (or 
source  elements),  which  are  suitable  assumptions  for  the 
hypothetical  flow  simulating  that  around  a  solid  body  in  a 
uniform  stream.  The  integral  (4-19)  then  becomes  time 


independent. 
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The  flux  is  also  expressed  by  the  velocity  on  S  as 


Q  =  Wcv*-  v)-'nd3 

$  w 


(4-20) 


Lagally's  theorem  for  the  force  on  singularities  gives 


and 


(4-21) 

Vs  v,  st 

~  (4"22) 

T,  Vs 


where  U«,  is  incoming  flow  and  Di  is  velocity  induced  by 
other  singularities  (vortices)  and  by  sources  if  any  at 
outside  S.  V~s  is  the  position  vector  of  the  source.  The 
second  term  on  the  right  of  each  equation  is  due  to  the 
motion  of  the  sources.  The  forces  which  act  between  the 
sources  inside  S  cancel  out  and  do  not  remain  in  the  total 
force. 

Substituting  Equations  (4-21)  and  (4-22)  into  (4-17) 
and  (4-18),  we  get 

+  ^ap  “  Pl*  =  " P \\(VU* ~  -  r%^dV  +  ?%ji(lX.+  vZ)*V 

Sw  v  v* 

+  (4-23) 

.  S  Vs 

and 

Sw  v 

% 

-f  $r5*C(LUt  Vi)f<*vs  (4-24) 

ys 


r 
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Construction  of  Hypothetical  Flow  Different  from  Betz's 
method  where  measured  pressure  automatically  gives  the 
hypothetical  flow,  it  is  necessary  to  construct  the 
hypothetical  velocity  in  the  present  method. 

As  mentioned  in  the  previous  section,  most  convenient 
is  that  the  hypothetical  flow  gives  the  same  lift  as  the 
real  flow.  This  requires  that  the  trailing  vorticity  be 
identical  in  both  flows.  Hence  there  is  a  possibility  that 
we  are  able  to  get  the  hypothetical  flow  by  measuring  the 
trailing  vorticity  and  using  Biot-Savart' s  law.  However,  it 


may  be  tedious 

and  time  consuming 

to  do  so. 

- 

The  author  has  come  up  with  an 

idea 

to 

use 

Navier -Stokes 

equation  provided 

that  the 

real 

flow 

is 

satisfying  it. 

The  equation  is 

written  in 

a 

form 

with 

vorticity  expressed  explicitly, 

(4-25) 

where  G  is  potential  of  gravity.  Here,  the  viscosity  term 
was  neglected,  because  it  is  presumed  to  be  small. 
This  is  consistent  with  the  assumption  stated  in 

Equations  (4-1)  and  (4-2)  .  This  approximation  makes  the 
equation  appear  to  be  identical  to  Euler's  Equation. 
However,  note  that  the  vorticity  has  difference  as  will  be 
mentioned  later. 

The  hypothetical  flow  is  supposed  to  be  inviscid  and  to 


satisfy  Euler's  equation. 
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=  &*  *  Qjl 0*  (4-26) 

Now,  the  real  flow  in  the  wake  is  considered  to  have 
two  types  of  vorticity.  One  is  boundary  layer  type 
vorticity,  which  originates  in  the  boundary  layers  on  the 
body  surface.  The  other  is  trailing  vorticity  which  can  be 
modeled  by  inviscid  flow. 

Write 

<UJ>  =  Oxbr*  (4-2  7) 

where  (UJr  =  vorticity  which  constructs  trailing  vortex 
system,  hence  contributes  to  circulation 
and  hence  to  the  lift, 

<UJq  »  vorticity  other  than  <uj y, 

originating  in  the  boundary  layers 
on  the  blade. 

The  vorticity  in  the  hypothetical  flow  is  just  trailing 
vorticity.  There  are  no  vorticity  which  originates  in  t^e 
boundary  layers. 

<UJ*  =  OJr  (4-28) 

Subtracting  Equation  (4-25)  from  Equation  (4-26),  we  obtain 

•Hr -it  -  ( v*-v)*<uDr-zj*uj)B  (4-29) 

This  is  now  a  differential  equation  to  give  £>*,  the  boundary 
condition  being  TL? ■  TD  on  the  edges  of  the  wake. 

It  is  not  advantageous  to  solve  this  in  an  ordinary  way 
for  three  unknowns  of  the  velocity  conponents.  By  making  an 
assumption,  things  can  be  greatly  simplified. 


\ 


That  is  : 


The  difference  of  the  two  velocities  is  parallel  to 
the  trailing  vortex, 
or  Tp*  -  TP  /  /  <uj>r  • 

This  assumption  is  essentially  the  same  as  Betz  made  in  his 
hypothetical  flow,  and  is  necessary  to  maintain  the  same 
lift. 

Equation  (4-29)  then  reduces  to 

~  <^8  (4-30) 

4.2  Applicat ion  to  2-D  Hydrofoil 

Formula  To  check  the  validity  of  the  method  described 
above,  experiments  were  conducted  first  with  2-D  hydrofoils 
so  that  the  results  could  be  compared  with  other  data. 


Figure  4-3  2-D  Hydrofoil  * 


For  a  steady  2-D  hydrofoil  or  aerofoil,  in  a  uniform 


flow  of  unbounded  fluid,  Equation  (4-23)  reduces  to 


+  Ffto  -  w*  *  -  P -  v*v:)d$  -  p^cnj* 

where  is  the  induced  velocity  due  to  vorticity 

distribution  (foil)  in  the  hypothetical  flow.  Since  the 
vorticity  distribution  in  the  hypothetical  flow  is  intended 
to  be  identical  to  that  which  contributes  to  lift  in  the 
real  flow,  contains  the  same  force  as  fRL  .  Besides,  it 

has  the  lift  due  to  the  induced  velocity  from  the  extra 
source  distribution.  But  this  lift  exactly  cancels  out  with 
the  last  term  on  right  side  as  internal  force,  even  if  there 
are  extra  sources,  provided  that  those  extra  sources  are 
inside  the  control  surface  S.  Therefore,  we  have 

For  unit  length  of  span  (or  for  2-D  foil),  the  surface 
integrals  reduce  to  line  integrals.  Also,  making  use  of  the 
relation  (4-19)  and  (4-20) ,  we  have 

F«d  =  -P$  'V*ZJ*)d5  -  ?V*>\(.d>-V)'n4S 

The  drag  acts  parallel  to  x-axis.  Taking  x-component , 
we  have 

Fee*  =  -  P -  pu 

■Svy  Sw 

*)  Note  that  the  coordinate  system  is  defined  differently 
from  the  propeller  case. 


If  we  take  x=const.  line  as  a  part  of  the  control 
s urface(  1  ine)  downstream  the  foil,  we  have 

'n  -  a  ,  *  v%  ,  zrZ  =  V  *  . 

Hence, 

/•Mwi  .j  /-$wj  + 

Fbox  *-f  \  -PU^  (Vv-Ujfy 

V  (4-31) 

The  drag  coefficient  Cd  is  given  by 


i:  PUoJc 


(4-32) 


Construction  of 


Flow  Equation  (4-30)  becomes 


(tix+HfllMvl  v/Ogj ) 

In  2-D  steady  flow,  there  should  be  no  trailing  or  shed 
vorticity  in  the  wake.  All  that  is  there  is  the  boundary 
layer  type  vorticity  in  the  viscous  wake.  Therefore,  there 
is  no  need  to  discriminate  one  type  from  the  other. 

Taking  y-component,  we  have 

CV*  )  J  =  ^‘[‘2’ )]  + 

Integrating  with  y, 

U^Z  +  vf  =  27**  +  2  C  ViCOg  otfr 


where  ^  is  one  of  the  edges  of  the  wake.  The  integration 
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should  be  only  in  the  wake  of  the  foil,  since  the  vorticity 
is  zero  outside  it. 

Transform  the  coordinates  from  x,y  to  e,n  in  which 
C-axis  is  parallel  to  the  streamline  at  the  center  of  the 
wake.  Since  there  is  no  relative  motion  between  the  two 
systems,  the  absolute  velocity  does  not  change  through  the 
transformation,  giving 

V?  ♦  Vy  =  2 +  Uy 


The  equation  then  becomes 

•J[ 

Uf  +  +  Vy  -t-  2 

■if  ^ 

The  equat  ion  still  has  two  unknowns  and  2Jy  .  To  solve 

this,  we  need  one  assumption.  The  assumption  stated  before 
should  be  applied. 

In  this  particular  case,  it  comes  as 

Uy  -  Uy 

This  is  explained  from  the  following  reasoning. 

The  trajectory  of  the  wake  near  the  trailing  edge  is  bent 
due  to  the  presence  of  lift  (bound  circulation) .  To  keep 
the  same  lift,  the  geometry  of  it  has  to  be  the  same. 
This  means  that  the  difference  in  the  velocity  between  the 
two  types  of  flow  is  permitted  only  in  ,  in  other 

if 

words,  Z JZj  »  Uy  has  to  be  maintained. 

The  equation  simplifies  to 


(4-33) 


Since  the  right  side  is  known  (given  from  the  experiment) , 
is  determined.  Through  the  inverse  transformation, 
and  Vy  are  now  given. 

To  reduce  possible  numerical  errors,  a  further  step 
should  be  taken. 

Since 

oj  -  2V* 

' 

the  integral  on  the  right  becomes 

2 

=  2  C  -  ( z//-  ) 

3-, 

and 

+  2.J  -  V*  +  V*  itivi  (4-34) 

Experimental  Results  Wake  surveys  were  conducted  in  the  MIT 
Water  Tunnel  using  two  hydrofoils  s 

A)  NACA  0012  Basic  Thickness  form 
and 

B)  NACA  66-309  a»0.8  with  modified  L.E.  *) 

Both  foils  have  the  dimensions  of  9" (chord)  x  20" (span)  and 
conpletely  span  the  water  tunnel.  The  experimental 
conditions  are  shown  in  Table  4-1  below. 


*  The  wake  survey  for  Foil  B  was  done  by  Moas[19]. 
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Table  4-1  Experimental  Conditions 


The  wake  surveys  were  done  by  traversing  the  laser 
along  midchord  plane  at  sevral  x=*const.  lines.  Figures  4-4 
through  4-6  show  measured  and  hypothetical  velocities 
( x-component)  and  vorticity  distribution  computed  from  the 
velocity  data.  Figure  4-7  show  the  values  of  the  profile 
drag  obtained  at  various  measurement  points  in  x  direction. 
Ideally,  Cd  value  should  be  independent  of  the  measurement 
point,  unless  large  dissipation  takes  place.  In  fact,  the 
results  shown  here  give  fairly  constant  Cd  along  x. 

For  the  NACA  0012  foil,  Cd  at  zero  degree  of  angle  of 
attack  is  somewhere  between  0.006  and  0.01  from  reference 
[20].  The  values  measured  by  the  LDV  are  considered  to  be 
fairly  reasonable. 

For  the  NACA  66-309  foil,  the  author  has  not  found  the 
curves  for  Cd.  However,  extrapolating  from  the  chart  for 
NACA  66-209,  the  angle  of  attack  4.5  degree  gives  C*=0.6. 
Then  the  corresponding  Cd  is  somewhere  between  0.008  and 
0.011,  to  which  are  conparable  the  values  obtained  by  the 
LDV. 

The  above  examples  have  well  demonstrated  the  validity 


of  the  present  method 
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4.3  3-D  Hydrofoil 

Wake  survey  was  also  done  with  a  3-D  (i.e.  finite 
span)  hydrofoil.  SNAME  Keel  Model  No. 3  was  attached  to  the 
side  window  of  the  water  tunnel  and  the  velocity 
distribution  behind  the  foil  was  measured.  The 
characteristics  of  the  foil  and  the  conditions  of  the 
experiment  are  listed  in  Table  4-2. 


Table  4-2  Dimensions  of  the  Model 

and  Experimental  Conditions 


Identification 

SNAME  Keel  Model  No. 3 

Span  (or  Height) 

7.320" 

Chord  length 

8.000” 

Sweep  angle 

40deg 

Section 

NACA  63-010 

Angle  of  attack 

4deg 

Reynolds  No.  (chord)  1  x  10' 

The  measurement  points  were  taken  at  z=*-5.0,  -1.0  and 
-0.5  inches  at  x*-0.4",  the  origin  of  the  coordinates  being 
at  the  trailing  edge  of  the  tip.  Along  the  line  parallel  to 
y-axis  the  velocity  was  measured.  (See  Figure  4-8  for  the 
difinition  of  the  coordinate  system. )  Also  the  velocity  was 
measured  across  the  tip  vortex  core  center,  traversing  the 
LDV  parallel  to  y-axis  and  also  parallel  to  z-axis. 


Il 
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Figure  4-8  3-D  Hydrofoil  and  the  Coordinate  System 


Figures  4-9  to  4-12  show  thus  obtained  velocity 
distributions.  At  z=-0.5"  (Fig. 4-9),  the  velocity  defect  is 
shallower  than  others.  This  is  because  this  measurement 
point  is  further  than  other  points  from  T.E.  due  to  the 
sweep  of  the  foil. 

It  is  interesting  to  see  that  the  velocity  distribution 
across  the  tip  vortex  core  shows  double  dents  (Fifure  4-12). 
This  also  appears  in  figure  4-13  which  corresponds  to  the 
traverse  parallel  to  z-axis.  These  two  figures  are 
revealing  that  a  part  of  the  boundary  layer  on  the  blade  is 
'sucked  in'  the  tip  vortex  and  swirl  around  the  vortex  core, 
which  was  also  observed  in  the  propeller  flow.  If  we  follow 
the  swirling  path  of  the  velocity  defect,  on  these  figures 
(4-12  and  4-13),  we  see  that  the  defect  becomes  shallower 
and  wider  as  it  turns. 


The  profile  drag  was  confuted  at  each  point  using  the 
measured  velocity  defect.*  The  spanwise  distribution  of  the 
drag  (coefficient)  is  shown  in  Figure  4-14.  At  the  mid-span 
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range,  Cd  is  constant.  Near  the  tip,  it  decreases  once  and 


then  in  the 

tip  vortex. 

it 

becomes  larger. 

The 

drag 

coefficient 

obtained  at 

the 

tip  vortex  may 

not 

be  so 

accurate,  because  the  drag  is  calculated  as  sectional  drag 
(in  other  words,  2 -dimensional ly)  there.  To  obtain  more 
accurate  value,  the  momentum  defect  has  to  be  taken  for  the 
whole  tip  vortex  area.  However,  there  should  be  no  doubt 
about  the  fact  that  tip  vortex  region  contains  considerable 
amount  of  momentum  loss  -  drag.  Even  at  zero  lift  (i.e. 
without  tip  vortex) ,  it  was  shown  that  the  sectional  drag 
decreases  toward  the  tip  and  then  increases.  This  'tip 
effect'  could  give  substantial  amount  of  difference  to  the 
drag  value  from  2-D  case  [21]. 

4.4  Application  to  Propeller 

Formula  To  manipulate  Equations  (4-23)  and  (4-24)  we  need 
to  be  more  specific  about  the  control  volume  and  the 
singularity  distribution. 

As  a  control  volume,  choose  a  circular  cylinder  of 
radius  Rc  whose  axis  (center  line)  lies  on  x-axis.  Name  the 
three  surfaces  which  surround  the  control  volume  as  S,  ,  S* 
and  S3  as  shown  in  Figure  4-15.  The  radius  Rc  is  to  be 
taken  large  enough  so  that  the  propeller  is  inside  this 


*)  The  formula  for  the  profile  drag  of  a  2-D  hydrofoil  can 
be  also  applied  to  a  finite  span  foil  if  we  assume 

i)  Induced  velocity  is  perpendicular  to  the  inflow 
ii)  Streamline  in  the  wake  is  parallel  to  x-y  plane. 


control  volume. 


For  the  case  of  constant  speed  rotation  of  the 
propeller  in  a  uniform  stream,  more  modifications  are 
possible  on  some  of  the  terms  in  Equations  (4-23)  and 
(4-24)  . 


If  we  write  t)  =  Tu'C^^j  &',t)  +  y££  (£© 


then 

—  JL -o/(<y  *r  f)  T')  I  ■+• 

at:  dT V 1  a©*aF 

=  M'/  P  3SL' 

a'C4=t  c  a©' 

From  blade  fixed  frame  of  reference,  the  flow  looks  steady. 

This  gives  -i^-/  =  O  • 
vc  /c»t 

Hence 


at  6trae' 


(4-35) 


The  second  term  on  right  side  of  Equation  (4-23)  becomes 

=  -  PS.  (ft  |f ' 


Kr«' 


«  r 


-  -  ?ci 


due  to  periodicity  if  no  blade  is  in  the  way  of  integration 
path.  Or,  if  the  blades  are  in  the  way,  it  becomes 


C  V'  fKg 
0  *0 


C  *  C 

0-c  w  >e,*AQt 

~  O  -  X>'l ,  -*0-0  *  TO'/.  -  O 

'0»O  I6*IK 

-  -  u/.  +  v'L 

/e-*o  i6-i 


=  o 


where  e,  designates  the  angular  location  of  the  surface  of 
one  of  the  blades,  *et  is  the  thickness  of  the  blade  in 
angle  in  e 1  direction  at  radius  r  and  K  is  the  number  of 
blades.  This  result  is  obtained  from  the  no-slip  condition 


on  the  blade  surface 

and  from 

periodicity. 

Hence  the 

integral 

is  zero 

in  any  case. 

The 

first 

case 

applies 

to 

and 

this 

integral 

becomes 

zero. 

/  9U« 

is 

zero  because 

the 

inflow 

is  uniform. ) 


lv£o rdrdx  =  O 

V*  r 


Similarly, 

P®0"*  -  P£2SSJ&f,ri< 

v  v 

-  -  pS2$$(nrxa/]*'  o  ydrcix  =  o 


and 


O 


Now,  let  us  discuss  the  fourth  term, 


on  the 
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right  of  Equation  (4-23) ,  The  integrand  takes  different 
forms  depending  on  the  surface  along  which  the  integration 
is  to  be  done,  since  the  normal  unit  vector  becomes 
different . 

On  St  ,  the  normal  unit  vector  is 

O 

TS  =  S 

on  Sx  ,  it  is 

O 

n  =-a 

and  on  S3  ,  it  becomes 

n  -  <sr  =  cose'J'+  si^e' Ik' 

To  see  if  the  integral  vanishes  (which  is  the  most 
convenient  for  us),  we  need  a  more  specific  expression  for 
(ps  .  The  velocity  potential  for  the  source  distribution 
which  corresponds  to  K-bladed  propeller  is  given  by 


£  err  j _ fleers-,  _ 

Ys  £,)))  41t  fox-*,)1 1  Y7+Yi-2TYs  Cos  (e'-Gj  - 


Ysde^drid%i 


(4-36) 


where  (xs  ,  r3  ,  es’  )  is  the  location  of  the  source  element  and 
<3(xs»r3,0al  )  is  the  intensity.  The  integral  is  over  the 
volume  of  the  source  for  one  blade.  Note  that  this  is  the 
potential  in  the  fixed  frame  of  reference,  xre,  although 
expressed  in  terms  of  xre 1  . 

Then, 


9^  99' 


-  rrs*s ince-  Gi-S&')&(.'Xs,Ys>es')  JA<Av 

VA 


If  we  expand  the  denominator,  we  have 
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_ 1 _ 

CC-X-TTs^V^Vl*  -  -vers,  COi  (©'-©J.'.  ■) )  a/* 


tCJMf^VyW  J3/- 


r  (  1  +-§-0.  cos(©/-©/-2g^)^-i|;a,‘cos2C©-^-^)  +  - 


where  CL  s  — —  ™*— ,  was  used  for  convenience. 

ComcoV 


And,  if  we  use  the  following  formula, 


*£'(?) *£(£) 


cos"9  = 


jC»o 


U=£<n 


we  obtain 


COSCn-iA)  9 


m*  sir' C e-e/ -  ->  % Cx v, n , fj.  ) 

(<.'x-‘Xs')l-*-yVv's7  -  am  cos  ce'-e,  -  }  ]3A 


^  :  even 


n  :  odd 


•  tsssffl  ic-f 


Then, 


a?  coV«-j£Ke'-e;-^,]  de,'dr}d*i 


•0  K 


(eV~  ^  4  ?|c*»  *  » 


X 


Si  *(©- ©; -  2g£.)J  c(  e/  d  n  ct  ys 


It  is  shown  that 

E  3in(6'-©s  -  )  =  O 


* 

£  cos (cn-asiK 


*-!§-  si*  ikTCe'-©; )  ;$  ic^i,i,-') 


a-^d  •>$■ 
©tJ'erwise 


(See  Appendix  A.) 


Hence,  we  obtain 


»  “*« jj^fc,- J*I”  «*]♦[  lr  ?■  C"->  -r *i"i«e-4'>) 

(  g:Suek  tbit  M-2fl  =  ^K«F  i  ) 

Together  with  the  normal  unit  vector,  the  integral  becomes  : 
On  S,  and  S2  , 

A\  =  tf  CiSMiK<e--«;)Jd%'ds«*«lyrf«ar 


Change  the  order  of  integration  and  execute  the  integral 


with  0‘  first.  We  have 


Therefore, 


On  S,  , 


ra7t 

\  si"<jKCe'-&')de'- o 

"'p 

$§£’«*'$  =o 


aX 

4^  (c°se'/'  +  sin©'  Ik'  )  rde'cjy 

y  o 

^1# % ft  [ft  sin^<'°" )](c*e’j*si»e’i>  )tfejdrsdxsrdedr 

By  orthogonality  of  trigonometric  functions,  we  get 

)  «(©'-©,' )co se'dS'  -  O 

<■** 

J  sWkc© '-©*')  si*  e'de'  =  o 


1 
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r 


V 


Hence,  we  have 

Sj 

The  term  which  corresponds  to  this  in  the  r  oment 
equation,  ^5  -ff'-flT* « )4$  ,  is  handled  similarly. 

Since 

Of  x  mi  »  *r£e/  =  ¥  fC-s.'M  ej'+  cose'/fe.' )  on  sL  or  sz, 
and  v  =  o(£0'  -  X  j' +  Co$e' If  )  on  S3  , 


the  integral  is  similar  to  fj  $&.nc/£ 

*»  9C 

as  was  used  to  handle  that  can  be  used. 
Then  it  is  obvious  that 


ffl 


-\<b* 

-)d$  *  O 


Equations  (4-23)  and  (4-24)  reduce  to 


The  same  argument 


El'Fw-R*  =  -  P^?r-V,c -v  iv, ) rsde/drsdKs 

<&w 


W«u+n^-Mt=  -  Pj|i r*(vis* 


-rfu*)A$  -  ?\\\%  vrs*(-vj -*  v, ) rsdd'dr,dxs 

~rc 


In  examining  the  components  of  lifting  force  in  the 
hypothetical  flow,  the  same  argument  as  that  used  in  2-D 
hydrofoil  can  be  used.  in  this  case  however,  there  is 
trailing  vortex  extending  to  downstream.  This  has  to  be 
taken  into  account. 

The  lift  Fu  is  decomposed  as 

F'  =  F*  ♦  *  F * 
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where  is  the  extra  lift  due  to  extra  sources.  This  term 
is  cancelled  by  the  force  on  the  extra  sources  due  to  bound 
vortices  inside  S.  The  force  equation  reduces  to 

ffc  =-? \\(vv«-v*u:u$  -  fjjfac -Vt  Vi ) 

% 

where  Vi  now  represents  induced  velocity  due  to  trailing 
vortex. 

For  the  case  of  moment,  the  extra  sources  give  rise  to 
extra  moment.  The  moment  in  the  hypothetical  flow  is 
decomposed  as 

This  time,  the  last  term  does  not  cancel  with  the  moment  on 
the  extra  sources,  but  produce  coupling  moments.  If  we 
segregate  those  moments  as  AMs,  the  moment  equation  becomes 

fjj  -  ?]\)%rs*(-Tj.+vt)r,<leidriA7s+aH^ 

where  U;  is  the  same  as  the  force  equation. 

Since  the  wake  is  divided  into  K-identical  wakes  from 
the  blades,  we  can  redefine  Sw  as  the  wake  on  S  for  one 
blade  only.  Then,  each  surface  integral  will  change  to  the 
sum  of  K  integrals. 

Define  Sw  s  a  part  of  control  surface  S  which  lies  in  the 
wake  of  one  blade 

Vs  s  volume  of  source  distribution  for  one  blade. 


We  can  choose  the  radius  Rc  of  the  cylinder  (control 
volume)  large  enough  so  that  Sw  lies  on  S4 .  If  this  is  the 
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O 

case,  we  have  '’HZ  =  -  £  and  U*  *  -V*. 
Hence, 

=  -SS(vZ-v*)*& 

n  A, 

And  we  obtain 

^o«>  =  -PK^Cwon-tfulT)^  “ 

5  m/ 


-  dvs 

Jtlg0  =  -  PKj{trx( vu,-iA£)4$  -  pic vj  4\  %rs  <tvs 

-  PkgV&^x2y*)dTs  +  *MLS 

V, 


(4-37) 


(4-38) 


The  axial  component  of  the  force  and  moment  will  be 
used  later  in  the  drag  calculation.  They  are  given  by 

fw*=  pk$( -  pkvifS^*-^0^ 

-  S\>, 


Sw 


-  PKffi?  Woof'S 

f1»,=  fKfi v(u.  14- v'ol)J*  -  PK((fr5gy,edvs 

IW  ^ 


(4-39) 


(4-40) 


Construction  of  Hypothetical  Flow  for  Propeller  The 
measurement  point  of  velocities  is  fixed  in  space  while  the 
propeller  turns.  If  looked  from  the  propeller  attached 
frame  of  reference,  the  measurement  point  draws  a  circle. 
The  direction  of  rotation  of  the  measurement  point  along 
this  circle  is  opposite  to  the  direction  of  propeller 
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rotation.  The  following  discussion  applies  to  the  arc  drawn 
by  the  measurement  point. 

If  we  take  e-component  of  Eq.uation  (4-30)  ,  we  have 


dV, 

at 


*  *■?  k  (i  i  *>'0  =  It1  *  vh(-k  i  vi‘ )  -  ( ) 


The  time  derivative  is  now  replaced  by  that  w.r.t.  e1,  if 


we  make  use  of  (4-35)  . 

dVo  _  o  dUo* 

at  —  ~  oC  ^  Q, 


1FET 


ae- 


And  also  partial  differentiation  w.r.t.  e  is  replaced  by 
that  w.r.t.  8‘,  giving 

u’,+  rSi  J 

Here,  the  relation  l/e  =  LL' +  was  also  used.  The 


equation  reduces  to 


r  ae't  2  (  +^r  +  H'  H  ~  Y^ae'OjH'  })  ^r^Tf) 


Integrating  with  respect  to  e 1  and  multiplying  by  2,  we 
obtain 


i£+vr  +  v;,  = 


-2C(u,co#),-^a;f^yci0,(4-4i) 

e' 


This  equation  still  seems  to  have  three  unknowns.  But  if  we 
use  the  assunption  stated  previously,  which  is 

D*  -  X)  / /«Jf  ,  we  only  need  to  solve  for  one  component  in 
the  transformed  coordinates. 

Let  us  define  a  coordinate  system  o'-gne/  (-axis  being 
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parallel  to  trailing  vorticity  *0/» .  n  and  ;  axes  are 
rather  arbitrary  as  long  as  the  three  axes  make  an 
orthogonal  system.  This  coordinate  system  should  be 
considered  locally,  or  at  each  point  of  interest,  since  the 
orientation  of  the  trailing  vorticity  may  vary  with  radius. 
The  transformation  is  done  by  first  shifting  the  origin  o' 
from  o  to  the  point  of  interest  and  then  rotating  the  axes 
so  that  5-axis  be  parallel  to  the  trailing  vorticity.  There 
is  no  relative  motion  between  o-xre'  and  o'-5nc  systems, 
which  gives 

=  Uf+Uj  .  (4-42) 

Hence  Equation  (4-41)  becomes 

Vj  V*  =  z((  U*o;B  Ur 

e* 

But  from  the  assumption,  we  have 

• 

Then  we  obtain 

Vf  *  and  U*  *  Uy  . 

Finally,  we  have 

v*  =  z/$z  -  zUv^ou^-  Uyco^)  rale'  (4-43) 

e' 

Once  Uf  is  obtained,  then  through  inverse  transformation, 
-  „  « 

1/*,  2/y  and  l/9,  can  be  obtained. 

Induced  "elocity  The  presence  of  trailing  vorticity  gives 
rise  tc  induced  velocity  which  changes  the  inflow  to  the 


blade 


To  fully  evaluate  the  forces  in  Equations  (4-39)  and 
(4-40),  the  induced  velocity  and  source  distribution  q  have 
to  be  obtained.  Since  the  induced  velocity  is  caused  by  the 
presence  of  trailing  vortices,  the  strength  of  them  has  to 
be  first  obtained.  This  information  is  then  fed  into  the 
Biot-Savart  law  to  calculate  the  velocity  on  the  blade. 

FPV-10  by  Kerwin[6]  and  PSFFPV  by  Greeley[15]  are  ones 
which  do  such  computations.  These  programs  were  designed 
for  the  confutation  of  field  point  velocity  due  to  a 
propeller.  They  are  also  capable  of  separating  each  element 
of  velocity,  like  velocity  due  to  trailing  vortex  only. 

The  source  distribution  is  more  difficult  to  get.  To 
know  it  exactly,  one  has  to  solve  a  boundary  value  problem, 
which  is  fairly  involved.  Since  it  is  not  the  main  theme  of 
this  paper  to  solve  for  source  distribution,  a  reasonable 
approximation  is  to  be  used  to  simplify  the  confutation. 
One  such  approximation  could  be  to  use  the  mean  values  of 
the  induced  velocities,  assuming  that  the  variation  of  them 
across  the  source  distribution  is  relatively  small.  The 
third  term  of  Eq.(4-39)  and  the  second  term  of  Eq.(4-40),  on 
right  sides,  are  then  approximated  by 

-pkuJS -  pk 

T»  T. 

Profile  Drag  Coefficient  The  most  convenient  way  for  the 
drag  coefficient  to  be  specified  would  be  to  have  it  given 
at  each  radius.  This  requires  strip  theory  assunftion,  in 

i 
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other  words,  flow  on  the  blade  is  2-dimensional  along  the 
radius.  While  it  is  said  that  the  flow  on  the  blade  is 
fairly  complicated  in  3-dimensional  way [14],  it  is  still 
considered  to  be  convenient  to  define  Cd  at  each  radius.  Of 
course  we  could  deiine  Cd  for  the  whole  blade.  However, 
this  will  not  give  much  detail  on  the  distribution  of  drag 
along  the  radius . 

Let  the  drag  be  DR£ (r ) ,  which  is,  according  to 
propeller  convention,  parallel  to  the  resultant  inflow. 
Then  profile  drag  coefficient  Cd  is  given  by 


Cd(r) 


DRiCr) 


(4-44) 


where  c  is  the  chord  length  at  each  radius. 

Assuming  the  drag  to  be  a  function  of  r  correspondingly 
requires  that  force  and  moment  be  given  as  functions  of  r. 
Suppose 

Fv*  -  <♦-«> 

r  Jec  :  section*! 

Mm.  '  (4-46) 

then 


Free  =  )<*5  ~  PKVa  Wv*~V*)<tS 

5^ 


+■  PK 


C  V?-Vr)dS 


(4-^7) 
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Ms«=  PKvjcy,v.-  v!i£)4s  +  ?KwClfr (4.48) 

iw  Suv 

The  surface  integral  in  was  replaced  by  line  integral. 
The  approximation  stated  on  the  induced  velocity  is  used 
here,  and  the  extra  moment  is  neglected.  It  has  to  be 

kept  in  mind  that  the  radius  of  measurement  does  not 
necessarily  correspond  to  the  radius  in  question  on  the 
blade.  Proper  correction  has  to  be  made  for  the  actual 
trajectory  of  the  viscous  wake. 

DRt(r)  is  given  by 

D«uOO  =  +  HsgCJCl  Co^] 

=  ~j*r[  E[«cCO  (Sfn/3co5*'+  CoS/flsiny  ) 

i-  (cos^cosf-  SinftsirtX  )  j 

where  =  (3;  -  (S 

and 

fi  :  hydrodynamic  pitch  angle  without  induced  velocity 
f. 3-  s  "  "  with  induced  velocity 

included 


.  - -  -■  ---- 


•rrc 


Figure  4-16  Expanded  View  of  Velocity  Diagram 


But, 


where 


Write 


Hence, 


coS/3 


yffi 


Si*/3 


Ta 

(  "V*1  +  V*&*  )  ^ 


1 


t>  1  v. 


TC*(^.y  J 


J  “Va/wd 


,  advance  coefficient. 


T/rc(-£)  -  T'  /  then 


Cos^S  - 


Sin/3  =. 


J' 


(i  a-  J> 1  )  '/*. 
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DniCr) 


J_f  r  r  7  c0?y 

Klhs~iO 


s<vi  y 

( i  + 


and  hence 


Mstef  CoS  y  j'sinY  | 


Ca(r) = 


_ Dr  i(r) _ 

-L-  pc  r*£V«  +  T'z) 


-2  cos  y _ 

f  cy^KCl+J'2)^ 


)' 


Substituting  Equations  (4-47)  and  (4-48),  we  have 

C’“(r'>  ‘  crWuZ'A  <  -UX^O 

+  i£,(  I  -7ta»s)jc/s 


O  rW)J(  Wr  U* )c/S  ] 

5w 


(4-49) 


where  is  given  by 


tv**  =  —  ~ - Uf*>  1 T* — 

(^--Uix)J'l+T>-  Vim7’ 


Results  The  velocity  and  vorticity  distribution  obtained 
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frora  the  experiment  and  shown  in  Chapter  3  were  used  to 

compute  profile  drag.  First,  it  was  tried  to  compute  the 

hypothetical  velocity  using  Eq.(4-43).  Finding  the  edges  of 

the  wake  from  vorticity  distribution,  integration  was  done 

* 

with  boundary  conditions  Vt  =  2J$  at  the  both  edges  of  the 
wake.  One  successful  (which  means  that  the  integral  started 
one  edge  and  arrived  at  the  other  edge)  example  is  shown  in 
Figure  4-17. 

This  hypothetical  flow  seems  wiggly.  From  the 
comparison  with  2-D  hypothetical  flow,  it  is  suspected  that 
this  is  due  to  the  noise  in  the  vorticity  which  was  used  in 
the  computation.  As  the  vorticity  is  given  by 
differentiation,  it  becomes  very  sensitive  to  the  error  in 
the  data.  Often  at  other  radii,  the  computation  of  the 
hypothetical  flow  failed  presumably  due  to  the  accumulation 
of  error  from  vorticity. 

To  improve  the  situation,  it  is  desired  to  further 
modify  the  integral  in  Eq.(4-43),  or  to  improve  the  accuracy 
of  the  vorticity  measurement. 

A  first  approximation  of  the  boundary  layer  theory 
tells  us  that  static  pressure  is  constant  across  the 
boundary  layer  or  wake.  For  a  moderately  loaded  propeller, 
this  assumption  is  considered  to  be  valid.  If  this  is  the 
case,  the  hypothetical  velocity  is  also  obtained  as  constant 
across  the  wake. 

To  see  this,-  a  few  points  along  a  streamline  of  viscous 
wake  were  picked  up  as  control  points  in  the  drag 
calculation.  At  far  downstream,  the  above  assumption  is 
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considered  to  be  more  realistic.  If  the  values  of  the  drag 
obtained  at  points  closer  to  the  blade  do  not  differ  so  much 
from  the1  one  obtained  at  the  furthest  point,  then  the 
validity  of  this  assumption  is  implied. 

Figure  4-18  shows  an  example  of  straight  line 
hypothetical  flow.  In  figure  4-19  are  plotted  the  values  of 
the  profile  drag  thus  obtained  at  three  different  positions 
along  the  viscous  wake  streamline  for  Prop. 4381.  The 
agreement  of  Cd  value  between  the  points  is  quite  good, 
proving  the  validity  of  the  assunption. 

Figure  4-20  illustrates  the  coefficient  of  profile  drag 
conputed  from  the  measured  velocity,  under  above  mentioned 
assumption,  at  various  radii  for  the  two  propellers. 

Toward  the  end  of  the  series  of  experiments,  the  author 
had  an  opportunity  to  test  Prop. 4381  with  smoother  surface 
condition  than  before.  This  was  due  to  re-surfacing. 
Hence,"  the  conditions  of  the  blade  surfaces  are 


Prop. 4381 

Relatively 

smooth,  marked 

by 

X 

M  II 

With  rough 

L. E . ,  marked  by 

A 

Prop. 4383 

Relatively 

smooth,  marked 

by 

o. 

The  two  different  surface  conditions  of  the  Prop. 4381  give 
considerable  difference  to  Cd. 

The  propeller  rpm  for  these  points  is  900rpra.  At  0.7R, 
however,  the  results  at  600rpm  and  1200rpm  are  also  added. 

The  radius  vs.  Cd  curves  show  fairly  complicated 
behavior.  At  mid-radius  (about  0.6R),  the  Cd  values  are  not 
so  much  different  from  what  they  would  be  in  2-D  flow.  But 
toward  the  outer  or  inner  radii,  they  become  lower.  And 
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this  '2-D  like'  region  is  narrower  at  Prop. 4383  than 
Prop. 4381,  which  may  imply  the  skew  effect.  And  at  further 
outer  or  inner  radii,  the  values  go  up  again.  This  is 
particularly  so  at  Prop. 4381,  but  not  so  distinct  at 
Prop. 4383 . 

This  trend  of  Cd  reduction  toward  the  tip  is  consistent 
with  what  we  observed  in  3-D  hydrofoil.  Since  the  aspect 
ratio  of  a  marine  propeller  is  much  smaller  than  airplane 
wing,  the  region  of  influence  by  the  tip  may  be  larger. 

Toward  the  propeller  hub,  a  similar  effect  could  occur, 

\ 

since  the  hub  corresponds  to  the  tip  on  the  other  side  of  a 
wing.  Besides,  the  flow  in  this  region  may  be  more 
complicated  due  to  the  presence  of  the  hub  and  interactions 
between  the  blades. 

The  low  Cd  at  outer  radii  (>£. 7R)  for  Prop. 4383  may  be 
partially  due  to  possibly  nontrivial  turbulence  level  of  the 
flow.  If  so,  the  present  method  will  give  an  error  to  Cd, 
since  that  is  neglecting  the  turbulence  components  which 
appear  as  Reynolds  stress  type  terms  in  the  Cd  calculation. 
From  the  analogy  to  a  swept  back  wing,  it  is  likely  that  the 
transition  from  laminar  to  turbulent  flow  is  taking  place 
there  in  the  thickened  boundary  layer  due  to  spanwise 
(radial)  flow[22]. 

Since  the  flow  is  quite  3-dimensional  at  the  both  ends 
of  the  blade,  it  may  not  be  adequate  to  define  sectional 
drag  there.  For  a  more  accurate  estimation  of  the  drag,  it 
is  desired  first  to  take  the  whole  drag  which  is  counted  in 
the  tip  (or  hub)  vortex  and  redistribute  that  over  the  (part 
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of  the)  radius. 

The  drag  coefficient  is  based  on  undisturbed  inflow  to 
the  blade.  But  near  the  hub,  it  could  be  possible  that  the 
inflow  is  already  retarded  due  to  the  presence  of  the 
propeller  shaft  boundary  layer.  If  this  is  the  case,  the 
actual  inflow  velocity  should  be  used.  This  will  increase 
the  Cd  values  near  the  hub. 
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5.  Conclusions 

The  following  conclusions  were  obtained. 

1)  LDV  has  proven  itself  to  be  a  handy  and  yet  powerful  tool 
in  the  analysis  of  the  flow  around  a  propeller. 

2)  The  wake  from  a  loaded  propeller  blade  is  categorized 
into  two,  trailing  vortex  wake  and  viscous  wake.  Trailing 
vortex  wake  is  created  by  the  tip  vortex,  whereas  the 
viscous  wake  comes  from  the  blade  boundary  layers.  Since 
the  origins  are  different,  they  have  to  be  considered  as 
separate.  Trailing  vortex  wake  appears  on  the  suction  side 
of  viscous  wake.  This  is  particularly  so  at  outer  radii. 
They  almost  line  up  at  inner  radii. 

3)  Relating  to  2),  the  vorticity  in  the  wake  is  also 
separated  into  two,  trailing  vorticity  and  boundary  layer 
type  vorticity.  By  finding  the  orientation  of  the  trailing 
vorticity  from  the  velocity  components,  new  coordinate 
system  o'-5nc  was  constructed  in  which  5-axis  is  parallel 
to  the  trailing  vorticity  at  the  locus.  This  method  works 
fairly  well  to  separate  the  two  vorticities.  5 -component  of 
the  vorticity  («j<)  appears  to  be  an  even  function  of  6  ' 
which  we  expect  it  to  be,  and  the  second  component  (o^)' 
shows  a  typical  vorticity  distribution  in  the  boundary 
layers.  It  was  also  found  that,  at  some  radii,  the  complete 
separation  of  the  vorticities  is  not  possible  only  by  the 


-82- 


coordinate  transformation.  This  is  due  to  the 
'contamination'  of  the  viscous  wake  vorticity  into  the 
trailing  vorticity. 

4)  The  Propeller  Field  Point  Velocity  Propgram  (FPV-10) 
predicts  velocities  fairly  well  within  the  limit  of 
potential  flow.  Also,  it  can  be  expected  that  the  modified 
version  of  it  (PSFFPV)  improves  the  ability  of  FPV-10. 

5)  A  formula  for  the  profile  drag  (viscous  sectional  drag) 
was  derived  in  terms  of  the  velocities  (real  and 
hypothetical) .  And  a  method  to  obtain  the  hypothetical 
velocity  was  suggested.  It  was  shown  that  the  formula 
worked  perfectly  with  2-D  hydrofoils.  For  propeller  case, 
however,  due  to  the  accuracy  problem  in  vorticity 
measurement,  the  hypothetical  velocity  sometimes  has  a 
problem.  It  is  desired  either  or  both  to  improve  the 
measurement  accuracy  of  the  vorticity  and  to  further  modify 
the  formula  for  the  hypothetical  velocity  so  that  it  can  use 
the  velocity  itself  rather  than  the  vorticity.  It  was  also 

shown  that  a  simple  assumption  for  the  hypothetical  velocity 

* 

(  ZSf  *  straight  across  the  wake)  could  be  used  for  the 
computation  of  the  profile  drag. 

6)  Radial  distribution  of  the  profile  drag  obtained  from  the 
velocity  data  showed  fairly  complicated  behavior.  At 
mid-radius  range,  the  value  is  cl  'ser  to  what  it  would  be  in 
2-D  flow.  But  toward  outer  r  inner  radii,  Cd  becomes 


-83- 


sub  st  anti  ally  low.  It  was  implied  that  considerable  amount 
of  momentum  defect  in  the  blade  boundary  layer  is  pulled 
into  the  tip  vortex,  thus  making  the  drag  appear  lower  in 
its  vicinity.  The  sarnie  trend  is  reported  in  aerofoil  flow 
even  at  zero  lift. 

7)  The  concept  of  profile  (or  sectional)  drag  may  lose  its 
validity  near  the  blade  tip  (or  the  propeller  hub)  since  the 
flow  there  is  fairly  3-dimensional  due  to  the  swirling  of 
the  tip  vortex  (or  hub  vortex).  It  is  therefore  desirable 
to  establish  a  way  to  estimate  the  total  momentum  and 
pressure  loss  (drag)  which  is  involved  in  the  tip  vortex 
(hub  vortex)  and  redistribute  that  over  the  (part  of  the) 
radius  (or  span)  to  get  correct  profile  drag  distribution. 

8)  In  this  paper,  the  turbulence  of  the  flow  (in  the  wake) 
was  assumed  to  be  small  and  was  neglected.  However,  for 
more  detailed  study  of  the  flow  around  a  propeller,  and  for 
more  accurate  estimation  of  the  profile  drag  especially  for 


a  skewed  propeller. 

it 

would 

be  desirable  to 

know 

the 

turbulence  level. 

To 

that 

end,  the  author 

feels 

it 

necessary  first  to  establish  a  method  to  segregate  the  noise 
in  the  Laser  signal  from  the  turbulence  component  of  the 
flow.  A  flow  visualization  technique  using  oil  film  on  the 
propeller  blade  would  supplement  the  LDV  measurement  to  get 
the  turbulence  level[23]. 
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APPENDIX  A 


i)  21  CosCCw-atfVtf- si*C0- 
4«i 


$  ^4  =  ^1  U  =  i,2,  •••  ) 
vU  yc  *  2 

©fcKev’wi&e 


[Proof] 


Cosifys-^Ccp  -  £j|&))sin(0  -  ) 

= -^  [si»>  ^|^))  -  s.*»C<M-u-OC*-agft>)] 


Bien, 


YL  Si*((n-3l+i.X<(> *=  6^n  zL  £ 
*=i  *~l 


£  _  i.(jwl*U(0-  ) 


Write  ^-ai^L  *  and  write  $  -  2L 

•-t 

Then, 


-i«-^ 


-  (  e 


+  e-‘~^ 


.  arc* 
e— -V  ) 


*jck  -  arc 


/•  a '*•’ If '‘'"Te"  \ 

-  (  e  K  +  •••  +e  *  +e,  *  * ) 

-u*2&  -tm-T?-  -IT  mu 

=  e  v  -  e  e 


c  i-  &im^' )£  =  eim^  c  i -  e*™ ) 


Since  Q  *=i  (n^integer),  the  right  hand  side  is  zero. 


Hence,  «0  if  ^  ^integer,  j  say. 
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If  =  j ,  then 


S  =  Zs 


-ZTX-ifc  < 

=  L  l 


K 


The  same  is  true  with  the  second  term,  ^L^siw(Cn-i£-i.X3>- . 

It  is  only  when  K=2  that  n-2j  +l=jK  and  n-2$  -l=j’K  hold 
simultaneously.  Otherwise,  we  have 


COs((n-2G)(fJ>- S\n(4>  -  ■SLfe) 


r  ±  K  S.-n<^K'<()) 


if  *\-29  =•  ^Kf  i.  c  *=1.2,3- •  ) 

a«d  rf  (C  *  Z 


ii)Z.SinC  4-*^)  =0 
[Proof] 

This  corresponds  to  the  case  when  n-2  =0  in  the  above 

(  i  ).  Since  jK±l*»0,  this  is  the  second  case.  Hence,  the 
sum  is  zero. 


Q*  E  •  D  • 


Table  2-1  Data  Acquisition  System 
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Table  2-2  Geometry  of  The  Propellers 


Number  of  Blades 

Expanded  Area  Ratio 

Section  Meanline 

Section  Thickness  Distribution 

Design  J 
Design  K 


5 

0.725 

NACA  a=0  .8 

NACA  66  with  NSRDC  modi¬ 
fied  nose  and  tail 
0.889 
0.213 


r/R 

tan 

c/D 

t/c 

0.2 

1.8256 

0.174 

0.2494 

0.3 

1.3094 

0.229 

0.1562 

0.4 

1.0075 

0.275 

0.1068 

0.5 

0.8034 

0.312 

0.0768 

0.6 

0.6483 

0.337 

0.0566 

0.7 

0.5300 

0.347 

0.0421 

0.8 

0.4390 

0.334 

0.0314 

0.9 

0.3681 

0.280 

0.0239 

Propeller 

4381(Skew=0  deg) 

r/R 

(deg) 

P/D 

f  h 

0.3 

0.0 

1.3448 

0.0368 

0.4 

0.0 

1.3580 

0.0348 

0.5 

0.0 

1.3361 

0.0307 

0.6 

0.0 

1.2797 

0.0245 

0.7 

0.0 

1.2099 

0.0191 

0.8 

0.0 

1.1366 

0.0148 

0.9 

0.0 

1.0660 

0.0123 

Propeller 

4382(Skew=36 

deg) 

r/R 

(deg) 

P/D 

f  /c 

0.3 

4.655 

1.4332 

0.0370 

0.4 

9.363 

1 .4117 

0.0344 

0.5 

13.948 

1.3613 

0.0305 

0.6 

18.378 

1.2854 

0.0247 

0.7 

22.747 

1.1999 

0.0199 

0.8 

27.145 

1.111' 

0.0161 

0.9 

31.575 

1.027 J 

0.0134 
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Table  2-2  (Continued) 


Propeller  449 7 (Warp-36  deg) 


r/R 

(deg) 

P/D 

f  /c 

0.3 

4.655 

1  .4332 

0.0370 

0.4 

9.363 

1.4117 

0.0344 

0.5 

13.948 

1.3613 

0.0305 

0.6 

18.378 

1.2854 

0.0247 

0.7 

22.747 

1.1999 

0.0199 

0.8 

27.145 

1.1117 

0.0161 

0.9 

31.575 

1.0270 

0.0134 

Propeller 

4383( Skew=72 

deg) 

r/F 

(deg) 

P/D 

f  /c 

0.3 

9.293 

1.5124 

0.0407 

0.4 

18.816 

1.4588 

0.0385 

0.5 

27.991 

1.3860 

0.0342 

0.6 

36.770 

1.2958 

0.0281 

0.7 

45.453 

1.1976 

0.0230 

0.8 

54.245 

1 .0959 

0.0189 

0.9 

63.102 

0.9955 

0.0159 

Propeller 

4498 ( Warp=72 

deg) 

r/R 

(deg) 

P/D 

f  /c 

0.3 

9.293 

1.5124 

0.0407 

0.4 

18.816 

1.4588 

0.0385 

0.5 

27.991 

1 .3860 

0.0342 

0.6 

36.770 

1.2958 

0.0281 

0.7 

45.453 

1.1976 

0.0230 

0.8 

54.245 

1.0959 

0.0189 

0.9 

63.102 

0.9955 

0.0159 
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Table  3-1  Distance  from  T.E.  along  Helical  Line 
to  Each  Measurement  Point 


r/R 

Distance/Chord 

Prop. 4381 

Prop. 43 8 3 

(X=-0.333R) 

(X=-0. 488R ) 

0.30 

0.462 

0.860 

0.40 

0.412 

0.688 

0.50 

0.428 

0.567 

0.60 

0.  502 

0.483 

0.70 

0.639 

0.427 

0.80 

0.880 

0.407 

0.875 

Advone*  Coefficient  J 


V*/VA  Vr/VA  -Vl/VA 


9‘ 


Figure  2—7  Result  of  Open  Uater  Test 

Prop.  DTNSRDC4498  Run  at  900rpm  MIT  MHL  Sept5. 1980 


Prop. 4381  at  Design  J  x=- . 33R,  r*.7R  Average  of  500 


I 


Vr/V 


Vr/V 


MKHMBti 


Vr/V 


02- 


Flgure  3— 8c  F I e I d  Point  Velocity  CTangential  component) 
Prop. 4381  at  x=_0.333R<  r=0.400R  at  Design  J 


Prop. 4381  at  x*-0.333R,  r*0. 500R  at  Design  J 


72.8 


Vr/VA 


0 


Vr/VA 


Vr/V 


-Vl/VA 


08- 


Vr/V A 


0- 


Figure  3— 14a  Field  Point  Velocity  C Axial  component) 
Prop. 4381  at  x=— 0.333R,  r-0.800R  at  Design  J 


F I gur •  3- 1  4b  Field  Po  int  Velocity  (Radiol  component) 
Prop.  4381  at  x*-0.333R,  r**0.800R  at  Design  J 


!i 

[i 

i 


i 


Vl/VA 


Angle  (deg.) 

Velocity  (Radial  component!) 
r=0 . 875R  at  Design  J 


2*.«  43.2  67.6  72.6 

Anal*  Cd*0.) 

Velocity  (Tangential  component) 
■0 . 87SR  at  Design  J 


Vr/VA 


L  3- 


Flpore  Field  Point  Velocity  CAxiol  component) 

P re*-  4381  at  x=-0.333Ry  r=0.9S0R  at  Design  J 


Prop. 4381  at  x=-0.333R,  r®0.950R  at  design  J 


Vl/VA  >1  TJ  -n  Vr/VA 


igurc  3- 1 7b  Field  Point  Velocity  C Radial  componanl) 
rop.4383  at  x=-0.488R,  r=0.300R  at  Design  J 


Prop. 4383  at  x=-0.488R,  r“0.300R  at  Design  J 


I  I  I  I  I  I  I  I  I  I  I  I  I  ••  I  I  I  I  I  I  I 

0.0  14.4  20.0  43.2  67.6  73 

Artgltt  (dt(.) 

Figure  3- 18b  Field  Point  Velocity  (Radial  component) 
Prop. 4383  at  x“-0.488R,  r“0.400R  at  Design  J 


mamnr 
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Figure  3- I  8c  Field  Point  Velocity  CT  anganl i a  I  component) 
Prop. 4383  at  x=-0 . 488R,  r=0.400R  at  Design  J 


Figure  3— 1 9a  Field  Point  Velocity  CAxial  component) 
Prop. 4383  at  x=-0.488R,  r=0.500R  at  Design  J 


L  9- 


F i gure  3- 20a  Field  Point  Velocity  (Axial  component) 
Prop. 4383  at  x=_0 . 488R„  r=0.600R  at  Design  J 


Figure  3-20b  Field  Point  Velocity  (Radial  component) 
Prop. 4383  at  x*~0.488R,  r*0.600R  at  Design  J 


Vr/VA 


Vr/V 


26.6  43.2 

tngl«  Cd*g  J 


FIjvjp*  3-22b  F I  •  I  d  Point  Velocity  CRodiol  componont) 
Prop. 4383  at  x=-0.488R,  r"0.700R  at  Dos i go  J 


-Vl/VA 


L  2  6- 


F  i  gore  3— 24c  Field  Point  Velocity  CTangentiol  component!) 
Prop. 4383  at  x=-0 . 488R,  r=0  875R  at  Design  J 


Figure  3—25  Vort icily  Distribution 

Prop. 4381  at  x=-0 . 333R,  r=0.700R  at  Design  J 


72.  e 


VA/ 


Ang I «  (d«g . ) 

Figure  3— 30a  Field  Point  Velocity  CAx i a  I  component ) 
Prop. 4381  at  x=-1  .000R,  r=0.700R  at  Design  J 


8.8  M.4  28.8  43.2  67.6  72.8  86.4  188.8  115.2  129.6  144.8 


Ang I •  (dag.) 

Figure  3— 30b  Field  Point  Velocity  CRodial  component) 
Prop.4381  at  x'=-1.000R,  r*0.700R  at  Design  J 


-0.48@~J 


-0.500—1 


-0 . 800— p-rrq~rrrr ;  tt  i  i  [  1 1 1  r  jTi  rr  ]~rrn  |'r  i  i  tjttttjtttt jtt  irj 

0.0  14.4  28  8  43.2  57.6  72.0  86.4  100.8  115.2  129.6  144.0 

Angl*  Cd«g . ) 

Figure  3—30C  Field  Point  Velocity  CTangential  component) 
Prop. 4381  at  x=-1.000R,  r-0. 700R  ot  Design  J 


Figure  3— 31a  Field  Point  Velocity  C Axial  component) 
Prop. 4381  at  x=-2.000Ry  r“0.700R  at  Design  J 


Vk/VA  Vr/VA  Vt/VA 


e. 


Figure  3-36  Field  Point 
Prop. 4381  at  x=-0 . 333R, 


Velocity  (Moosured) 
r=0.700R  at  Design  J 


Figure  3-37  FPV-10  Output  ot  x— 0.470R,  r-0.924R  Prop. 43 83 
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Figure  3—40  Transformed  Vort icily 

Prop.4381  al  x=-0.333R<  r=0  300R  ol  Design  J 


Figure  3-41  Transformed  Vorlieily 

Prop.4381  al  x*-0.333R,  r*0.400R  al  Design  J 


72.  • 


Vert  telly  *>/<VA/R)  *0  "H  Vorl  let  ly  *»/<VA/R> 
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igure  3—42  Transformed  Vorticity 

rop.438t  at  x=-0.333Ry  r=0.500R  at  Design  J 


Figure  3-43  Transformed  Vorticity 
Prop. 4381  at  x=-0.333R,  P-0.600R 


at  Design  J 


at  Design  J 


at  Design  J 


Vorl  lolly  «*V<VA/RJ  -q  Vorllelly  "/tVA/RJ 


-140- 


igur«  3—48  Transformed  Vort icily 

rop.438t  at  x=-0.333R<  r=0.87SR  at  Design  J 


Figure  3—49  Transformed  Vort Icily 
Prop. 4383  al  x*-0.488R,  r*0.300R 


al  Design  J 


Vartlel 


Anal*  Cdag.) 

Figure  3—50  Transforuad  Vort icily 

Prop. 4383  at  x=-0.488R,  r®0 . 400R  at  Design  J 


Figure  3-51  Transformed  Vorlicily 
Prop . 4383  ol  x«-0.488R,  r-0.500R 


ol  Design  J 


a l  Design  J 
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